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Long-Range  Surface- Plasma  Wares  on  Very  Thin  Metal  Films 

Dror  Sarid 

Optical  Sciences  Cantor.  University  of  Arizona,  Tucson,  Arizona  85721 
(Received  17  August  1981) 

The  dispersion  aqusttoa  of  Injected  surface-plasms  waves  that  propagate  on  thin 
metal  Sims  has  been  solved  as  a  function  of  the  film  thickness,  and  splitting  of  the 
modes  Into  two  branches  Is  observed.  For  one  branch  the  Imaginary  part  of  the  prop¬ 
agation  constant  goes  to  zero  as  the  thickness  of  the  metal  decreases.  Reflectivity  cal¬ 
culations  agree  with  this  result,  which  predicts  that  one  can  obtain  propagation  dis¬ 
tances  that  are  more  than  1  order  of  magnitude  larger  than  observed  before. 

PACS  numbers :  73.S0.Dt 


The  lifetime  or  decay  constant  of  a  surface- 
plasma  wave  (SPW)  that  propagates  on  a  metal 
surface  is  an  important  physical  quantity  because 
it  is  a  strong  function  of  the  properties  of  the 
metal  surface,  and  can  therefore  serve  as  a  sen¬ 
sitive  probe  for  characterizing  surface  conditions 
such  as  roughness  and  composition.1*4  There  is 
only  one  SPW  mode  for  thick  metal  films,  and 
the  upper  limit  of  its  propagation  range  is  de¬ 
termined  by  the  complex  refractive  index  of  the 
metal  and  the  loss  introduced  by  the  medium 
bounding  the  metal.’ 

The  existence  of  two  thermally  excited  SPW 
modes,  one  symmetric  and  one  antisymmetric, 
that  propagate  on  an  unsupported  thin  metal  film 
(the  two  media  bounding  the  metal  film  are  identi¬ 
cal)  has  been  discussed  theoretically  and  verified 


experimentally  by  using  electron  scattering  tech¬ 
niques.8’7  The  theory  that  predicted  the  existence 
of  these  two  modes  did  not  treat  the  problem  of 
their  lifetime  and  optical -wavelength  dependence, 
and  since  the  resolution  and  signal-to-noise  ratio 
of  the  experiments  were  low,  the  linewidth  of 
these  modes  could  not  be  measured. 

A  comprehensive  analysis  of  the  dispersion  of 
SPW  that  propagate  on  various  combinations  of 
thin  metal  films  sandwiched  between  thin  dielec¬ 
tric  films  revealed  the  existence  and  splitting 
of  the  SPW  rfiodes  as  the  metal  thickness  de¬ 
creases.*’9  The  theory  dealt  only  with  the  depen¬ 
dence  of  the  real  wave  vector  on  the  real  part  of 
the  frequency,  and  did  not  treat  the  properties  of 
the  Imaginary  part  of  the  frequency,  which  is 
associated  with  the  lifetime  of  the  thermal  SPW. 
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The  theory  of  the  dispersion  and  lifetime  of 
thermal  SPW  that  propagate  on  supported  or  un¬ 
supported  thin  metal  films  has  been  discussed  by 
Fufcui,  So,  and  Normandin.10  Their  theoretical 
conclusion  was  that,  for  an  unsupported  metal 
film,  there  is  a  SPW  mode  that  has  a  lifetime 
which  increases  as  the  thickness  of  the  metal  film 
decreases,  while  for  the  supported  metal  film 
such  a  mode  does  not  exist. 

In  this  Letter  a  different  case  is  treated,  name¬ 
ly  that  of  “Injected"  SPW,  where  the  frequency  is 
a  real  quantity,  while  the  wave  vector  is  complex. 
It  is  shown  that  for  both  the  supported  and  unsup¬ 
ported  metal  films  one  can  excite  a  SPW  mode 
having  a  decay  constant  that  goes  to  zero  as  the 
film  thickness  becomes  small  enough,  and  that 
this  mode  can  exist  even  in  the  visible  range  of 
the  spectrum.  The  reflectivity  at  the  base  of  a 
prism  that  is  placed  close  to  a  thin  metal  film 
has  also  been  calculated  with  identical  results. 
The  theory  of  the  splitting  of  the  SPW  modes  and 
their  properties  and  means  of  excitation  is  dis¬ 
cussed;  this  long-range  mode  offers  a  unique  op¬ 
portunity  for  probing  the  characteristics  of  metal 
surfaces. 

The  geometry  of  the  thin  metal  film  that  is 
coupled  to  a  prism,  and  that  supports  the  SPW, 
is  shown  in  Fig.  1,  where  nw  j*t  + 
and  are  the  refractive  indices  of  the  substrate, 
metal  film,  gap,  and  prism,  respectively,  and  t 
and  s  are  thicknesses  of  the  metal  and  the  gap, 
respectively. 

The  standard  dispersion  relation  for  a  free  TM 
mode  {s»\)  that  is  guided  by  a  thin  dielectric, 
semiconductor,  or  metal  film,  is  given  by11 

KBt  »tan"  l(klQ)  +tan*  ,  (1) 

where  K  -  2v/\  X  is  the  free-space  optical  wave- 


FTr?.  1.  The  geometry  of  the  prism  to  surface- 
plasmoa  wave  coupler.  Note  the  antisymmetric  (trans¬ 
verse)  electric  held  distribution  of  the  long-range 
mode. 
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length,  *,0  *  tA/nJB,  and  klt  •  Here  A* 

■  ,*  -0*,  C,»0*-ii,*,  A3  is  the  prop¬ 

agation  constant  of  the  guided  wave,  <■«,',  and 
m  is  the  mode  number.  Equation  (1)  is  an  implic¬ 
it  complex  equation  in  the  complex  variable  0 
*  0'*  i  0*.  where  3'  is  related  to  dm  speed  of  the 
mode  by  v*  c/K0‘ ,  and  3*  is  related  to  the  ampli¬ 
tude  attenuation  of  the  mode  by  a  »A0*. 

A  considerable  simplification  In  the  solution  of 
Eq.  (1)  can  be  made  by  defining 

*  »(*l0  - 1  X*„  -  <)/(*„ .  i X*u  * « ).  ( 2) 

From  Eqs.  (1)  and  (2)  one  obtains  the  real  implic¬ 
it  dispersion  relation  with  real  variables 

,  tan((3'/B-)lnJ*|  ]♦*-/*'  -0,  (3) 

where  k'  and  k“  are  the  real  and  imaginary  com¬ 
ponents  of  k,  respectively.  To  solve  Eq.  (3),  one 
chooses  a  value  for  0',  and  then  searches  tor  the 
value  of  0*  such  chat  Eq.  (3)  is  satisfied.  Conse¬ 
quently,  one  has  to  use  complex  numbers  only  for 
the  evaluation  of  A,  B,  C,  and  A;  a  simple  com¬ 
puter  routine  will  solve  for  0#.  The  film  thick¬ 
ness  is  obtained  by 

ln|Jr|  /2KJT.  (4) 

All  the  modes  we  are  considering  here  have  a 
mode  number  m  =  0.  By  using  Eq.  (1)  it  is  possi¬ 
ble  to  generate  curves  of  the  Interdependence  of 
the  three  variables  0',  0*  and  t,  for  a  given  set 
of  refractive  indices.  These  calculations  have 
been  carried  out  for  various  combinations  of  pa¬ 
rameters  at  optical  wavelengths  extending  from 
the  visible  to  the  far  infrared. 

For  a  thick,  unsupported  metal  film,  in  which 
the  two  bounding  media  are  identical,  the  speed 
and  decay  constant  of  the  SPW*s  on  each  one  of 
the  interfaces  are  the  same.  Each  one  of  these 
two  SPW' s  has  an  evanescent  wave  extending  both 
into  the  dielectric  medium  and  into  the  metal 
film.  If  the  metal  film  is  thick  enough,  the  eva¬ 
nescent  waves  inside  the  metal  that  belong  to  the 
two  SPW's  do  not  overlap:  As  the  film  thickness 
decreases,  the  evanescent  waves  of  the  otherwise 
decoupled  modes  begin  to  overlap,  and  a  trans¬ 
verse  standing  wave  is  established.  The  degen¬ 
erate  SPW  mode  therefore  splits  into  one  sym¬ 
metrical  and  one  antisymmetrical  mode  (referring 
to  the  transverse  electric  field  distributions). 
These  symmetric  and  antisymmetric  modes  are 
characterized  by  a  range  that  increases  and  de¬ 
creases,  respectively,  as  the  film  thickness  de¬ 
creases. 

For  a  thick,  supported  metal  film,  in  which  the 
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two  bounding  modi*  are  not  identical,  the  modes 
that  propagate  on  each  one  of  the  two  interfaces 
have  different  speeds  and  decay  constants  and 
their  evanescent  fields  inside  the  metal  do  not 
overlap.  As  the  metal  thickness  decreases,  these 
evanescent  fields  begin  to  overlap,  and  each  one 
of  these  two  modes  splits  into  two  branches,  as 
in  the  unsupported  case,  one  having  a  decreasing 
range  and  one  having  an  increasing  one.  In  the 
supported  metal  film  case,  however,  extra  inter¬ 
fering  terms  between  the  two  evanescent  waves 
inside  the  metal  will  limit  the  range  of  the  anti¬ 
symmetric  mode  if  the  difference  between  the  re¬ 
fractive  indices  of  the  two  media  bounding  the 
metal  is  too  large. 

In  Fig.  2  the  dependence  of  3'  (dashed  lines)  and 
3#  (solid  lines)  on  the  metal  thickness  t  is  shown 
for  the  case  of  the  free  SPW.  The  parameters 
used  are  *„■  1.5,  ■0.0657*  *4  (Ag),  *,*1.55, 

and  A  *6328  A.  One  observes  the  decrease  in 
both  3'  and  3*  for  the  antisymmetric  mode  as  t 
decreases.  The  physical  validity  of  the  results  of 
these  calculations  has  been  verified  by  checking 
the  value  of  the  evanescent  propagation  constants, 
A  and  C,  for  each  resultant  3',  3",  and  f,  to  in¬ 
sure  that  the  SPW  is  actually  confined  to  the  thin 
metal  film.  For  each  case  the  magnitude  and 
phase  of  the  electric  and  magnetic  fields  across 
the  thin  metal  film  has  been  calculated.  The 
absolute  value  of  the  (unsealed)  magnetic  field 
Hy(x)  is  shown  in  Fig.  3  for  the  metal  film  that  is 


0  500  1000 


t(S) 

FIG.  2.  The  real  (KB')  and  imaginary  (KB")  propaga¬ 
tion  constants  of  a  surface-plasma  wave  propagating 
on  a  thin  Ag  film.  Here  a,*  1.5,  *3  *  1.55,  A ■  6028  k, 
and  s »  A.  Note  that  the  decay  constant  a*  KQ"  goes  to 
zero  as  the  metal  film  thickness  decreases. 


in  the  y-a  plane  with  a  thickness  f  *  200  A.  In  (a), 
3*1.5507  +  10.0001615,  A* 0.39  +  1 0.00064,  and 
C  *0. 047  +  i  0.0053.  Outside  the  film,  the  field  is 
evanescent,  while  in  the  film  it  obtains  a  small 
dip.  In  (b),  3*2.14  +  10.0325,  A*  1.53  +  10.046, 
and  C*  1.44+  10.059,  and  the  field  obtains  a  value 
of  zero  at  about  the  center  of  the  film.  There  are 
two  electric  fields,  Ea{x)  and  £,(x),  the  first 
being  proportional  to  Hr{x)  in  each  one  of  the 
three  regions  and  for  the  two  modes.  The  second 
field  Is  proportional  to  the  derivative  of  B,(x) 
with  respect  to  x  in  each  one  of  the  three  regions, 
with  (a)  and  (b)  interchanged.  One  observes  that 
In  (a)  the  evanescent  fields  penetrate  the  media 
bounding  the  film  much  deeper  than  in  (b).  It  is 
assumed  that  the  only  loss  the  SPW  suffers  Is 
due  to  the  dissipation  of  power  inside  the  metal 
film.  Since  the  mode  in  (a)  has  a  smaller  fraction 
of  its  field  inside  the  metal  than  mode  (b)  has,  it 
will  have  a  larger  propagation  range. 

In  order  to  assess  the  feasibility  of  exciting 
this  long-range  mode,  the  reflectivity  at  the  base 
of  the  prism,  as  shown  in  Fig.  1,  has  been  cal¬ 
culated  for  the  same  material  parameters  used 
for  the  free  SPW.  It  was  found  that  in  order  to 
observe  the  extremely  narrow  dip  in  the  reflec¬ 
tivity,  which  is  associated  with  the  excitation  of 
the  loug-range  mode,  one  has  to  scan  across  the 
angle  of  incidence  of  the  input  beam  with  a  very 
high  resolution.  Since  this  mode  has  a  large 
evanescent  field,  the  angle  for  optimum  excita¬ 
tion  is  close  to  the  critical  angle.  As  a  result, 
small  changes  in  the  incidence  angle  produce 


FIG.  3.  The  abeolute  value  of  the  transverse  mag¬ 
netic  field  distribution  of  the  (a)  long-range  and 
(b)  short-range  modes  (not  to  scale.) 
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large  changes  in  the  penetration  depth  of  the 
evanescent  fields.  For  each  angle  of  incidence, 
therefore,  die  gap  thickness  s  has  to  be  adjusted 
so  that  the  proximity  of  the  prism  will  have  ex¬ 
actly  the  same  damping  effect  on  the  otherwise 
freely  propagating  SPW.  If  s  is.  too  large,  the 
excitation  will  be  inefficient,  while  if  s  is  too 
small,  the  proximity  of  the  prism  will  load  and 
wipe  out  the  SPW. 

Reflectivity  calculations  with  use  of  the  same 
parameters  as  in  Fig.  2  yield,  for  a  metal  thick¬ 
ness  of  200  A,  the  following  results:  The  critical 
angle  is  32.6244°,  the  angle  at  maximum  absorp¬ 
tion  is  32.6414°  and  the  half -width  angle  is  0.004°. 
The  distance  which  the  long-range  mode  propa¬ 
gates  on  the  200  A  thin  film  until  its  power  de¬ 
cays  by  a  factor  of  1/e2  is  therefore  300  Mm, 
which  is  27  times  larger  than  for  a  SPW  that  prop¬ 
agates  on  the  surface  of  a  thick  Ag  slab.  Since 
the  range  of  this  mode  is  very  large,  its  study  in 
the  case  that  the  surfaces  are  randomly  rough 
can  provide  useful  quantitative  information  about 
roughness-induced  attenuation.  As  a  result,  this 
mode  can  be  used  as  a  probe  for  characterizing 
the  properties  of  surfaces  and  thin  metal  films. 
This  long-range  mode  can  also  be  utilized  for 
nonlinear  interactions  such  as  second-harmonic 
generation  where  the' nonlinearity  is  introduced 
by  the  metal  itself  or  by  the  bounding  media  via 
the  evanescent  fields.  By  imbedding  the  thin  met¬ 
al  film  in  an  (optical)  power-dependent  refractive- 
index  semiconductor,  one  obtains  for  the  long- 
range  SPW  a  power -dependent  propagation-con¬ 
stant  which  can  be  used  for  bistability  studies.1* 

In  summary,  this  theory  predicts  that  it  is  pos- 
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sibte  to  choose  material  and  geometries!  parame¬ 
ters  such  that  a  long-range  SPW  mode  could  be 
excited  on  a  thin  metal  film,  the  range  being 
more  than  1  order  of  magnitude  greater  than  ob¬ 
served  before.  This  mode  can  be  -used  for  the 
study  of  rough  surfaces  and  various  nonlinear  in¬ 
teractions  where  a  large  interaction  range  is 
desirable. 
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Consider  a  TM  polarized  monochromatic  wave  that  is  in¬ 
cident  on  the  base  of  a  prism  from  which  it  is  totally  internally 
reflected.  If  a  metallic  substrate  is  placed  in  the  vicinity  of 
the  prism  so  that  the  evanescent  Helds  of  the  reflected  light 
penetrate  the  metal,  and  the  incident  angle  is  appropriate,  a 
surface-plasma  wave  can  be  excited  along  the  metal  gap. 
Otto1  has  analyzed  this  geometry  and  has  shown  both  theo¬ 
retically  and  experimentally  that  by  measuring  the  reflectivity 
of  a  plane  wave  as  a  function  of  the  angle  of  incidence  one 
obtains  both  the  speed  and  the  range  of  the  surface-plasma 
wave  from  the  position  of  the  minimum  and  the  width  of  the 
Loren  tzian  absorption  curve,  respectively.  In  his  geometry 
the  metal  thickness  was  much  larger  than  the  decay  length  of 
the  surface- plasma  field  in  the  metal.  Kretschmann2  con¬ 
sidered  a  different  geometry  in  which  a  thin  metal  film  is  de¬ 
posited  on  the  bottom  of  the  prism  and  a  plane  wave  incident 
on  the  base  of  the  prism  excites  a  surface-plasma  wave  on  the 
lower  side  of  the  metal  film.  The  excitation  in  this  case  takes 
place  via  the  optical  fields  that  penetrate  the  metal.  As  with 
the  Otto  geometry,  the  excitation  is  possible  because  the  speed 
of  the  surface-plasma  wave  is  slower  than  the  speed  of  the  light 
in  the  prism. 

Weber  and  Ford3  have  recently  shown  that  the  fields  ac¬ 
companying  the  surface-plasma  wave  are  larger  than  those 
associated  with  total  internal  reflection  in  the  absence  of 
surface- plasma  waves.  In  their  theory  they  assume  a  general 
coupling  mechanism  that  is  100%  efficient  and  equate  the 
optical  power  density  injected  by  the  coupler  under  steady 
state  conditions  to  the  power  dissipated  by  the  surface-plasma 
wave. 

They  show  that  the  excitation  of  a  surface-plasma  wave  at 
an  optical  frequency  can  give  rise  to  an  enhancement  of  the 
optical  electric  Helds  at  a  silver  metal  surface  by  a  factor  as 
large  as  60.  Their  theory,  however,  is  approximate  because 
in  calculating  the  enhanced  Helds  they  use  the  free  surface- 
plasma  Helds  instead  of  the  true  Helds  in  the  presence  of  the 
excitation  medium  (a  prism  or  a  grating,  for  example.) 

In  this  Letter  we  wish  to  show  theoretically  that  (1)  it  is 
possible  to  obtain  an  exact  solution  to  the  Held  enhancement 
for  both  the  Otto  and  the  Kretschmann  geometries  and  (2)  our 
new  geometry4  makes  it  possible  to  obtain  Held  enhancements 
that  are  as  much  as  an  order  of  magnitude  larger  than  those 
obtainable  with  the  two  other  geometries. 

The  three  prism -coupling  geometries  are  shown  in  Fig.  1, 
where  I  and  III  correspond  to  the  Otto  and  Kretschmann 
geometries,  respectively,  and  II  to  our  geometry.  In  each  case 
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Fig.  1.  Three  prism-coupling  geometries,  for  exciting  a  surface- 
plasms  wave:  I,  IH,  and  Q  correspond  to  the  Otto,  Kretschmann,  and 
our  geometries,  respectively. 

t  denotes  the  metal  thickness,  and  the  black  dot  in  each  di¬ 
agram  indicates  the  point  where  the  optical  Helds  are  evalu¬ 
ated.  In  geometries  I  and  II  the  evanescent  Helds  created  at 
the  prism  surface  need  to  penetrate  an  intermediate  medium 
separating  the  prism  and  the  metal  film.  The  exponential 
decay  of  the  Held  amplitude  in  this  medium  is  expressed  by 
the  factor  h  *  exp(-KCS),  where  K  *  2r/X,  X  is  the  optical 
wavelength  in  free  space,  S  and  n<i  are  the  thickness  and  re¬ 
fractive  index  of  the  gap,  C2  *  d2  -  r»i  and  K&  is  the  propa¬ 
gation  constant  of  the  surface- plasma  wave  which  has  only  a 
real  part  for  the  plane  wave  approximation.  In  this  work  we 
assume  a  varying  thickness  S  so  that  h  *  0.1  for  any  given 
angle.  The  reason  for  this  choice  is  that  under  these  condi¬ 
tions  the  reflectivity  obtains  its  minimum  value  for  the  ap¬ 
propriate  angle  for  most  of  the  cases  we  consider.  It  can  be 
seen  that  geometry  II  merges  into  the  Otto  geometry  for  t 
sufficiently  large  and  h  <  1  and  into  the  Kretschmann  ge¬ 
ometry  for  t  sufficiently  small  and  h  *  1.  For  geometry  II  to 
be  maximally  effective,  the  refractive  indices  of  the  media 
bounding  the  metal  Him  on  either  side  need  to  be  approxi¬ 
mately  the  same. 

Now  in  the  general  case  of  a  stratified  medium  characterized 
by  complex  dielectric  constants  the  reflectivity  of  a  TM 
polarized  monochromatic  plane  wave  at  the  surface  of  the 
medium  and  the  Helds  excited  in  each  of  the  media  can  be 
computed  by  standard  techniques.5  For  the  plane  wave  as¬ 
sumption  to  be  applicable  to  our  geometry  it  is  necessary  that 
the  range  of  the  surface-plasma  wave  be  appreciably  smaller 
than  the  width  of  the  incident  optical  beam.  In  the  visible 
range  of  the  spectrum,  where  the  latter  range  is  <1  mm,  this 
requirement  is  easily  satisfied,  and  the  simple  plane  wave 
calculations  yield  accurate  results  for  all  three  geometries. 

To  demonstrate  the  Held  enhancement  obtainable  with  our 
geometry  we  choose  parameters  appropriate  to  the  case  of  a 
silver  metal  film  having  a  dielectric  constant  given  by  - 16  4- 
i'0.52  at  a  wavelength  of  632.8  nm,  bounded  by  two  dielectric 
media  having  a  refractive  index  of  1.5.  Our  prism  is  assumed 
to  be  rutile  with  a  refractive  index  of  2.875.  Figure  2  show's 
the  amplitude  reflectivity  at  the  base  of  the  prism  for  three 
metal  film  thicknesses:  t  *  2000, 70,  and  50  nm  with  h  *0.1. 
Note  how  the  single  Lorentzian  absorption  curve  associated 
with  the  thick  metal  film  splits  into  two  Lorentzians  as  the 
film  thickness  decreases.  The  splitting  is  a  direct  result  of  the 
lifting  of  the  degeneracy  in  the  propagation  constants  of  the 
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Fig.  2.  Reflectivity  aa  a  function  of  the  incident  angle  inaide  the 

prism  for  our  geometry  where  n 0  »  n-i.  ( — >,  (---),  and  ( - )  are 

for  t  »  2000. 70,  and  50  nm,  respectively.  Note  the  splitting  of  the 
single  Lorentzian  into  two  Lorentzians  as  the  metal  film  thickness 
decreases.  The  wide  and  the  narrow  Lorentzians  correspond  to  the 
short-range  and  long-range  surface-plasma  waves,  respectively. 
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Fig.  3.  Dispersion  curves  of  the  free  surface-plasma  waves:  Im(j3) 
(dashed  line)  and  t  (continuous  line)  as  a  function  of  Re(/3).  Note  the 
divergence  of  £  for  Re(fl)  that  corresponds  to  the  thick  metal  case,  and 
the  two  values  of  Re(fl)  and  Im(/3)  for  each  t  for  the  thinner  metal 

case. 


two  independent  surface-plasma  waves  propagating  on  each 
surface  of  the  thick  metal.  Once  the  metal  film  thickness 
becomes  small  enough,  the  fields  of  these  two  surface-plasma 
waves  inside  the  metal  overlap,  and  the  degeneracy  is  lifted. 
For  such  a  thickness,  one  long-range  and  one  short-range 
surface- plasma  wave  can  be  excited  as  a  single  eigenmode  that 
satisfies  all  the  boundary  conditions  simultaneously.4  One 
observes  that  the  Lorentzian  width  at  half  maximum  for  the 
short-range  and  for  the  long-range  mode  increases  and  de¬ 
creases,  respectively.  The  width  consists  of  two  contributions, 
one  from  the  power  dissipation  in  the  metal  film  and  one  from 
the  prism  loading.  This  analysis  differs  from  that  of  Weber 
and  Ford,  because  it  takes  into  account  the  important  fact  that 
any  efficient  means  of  excitation  has  in  principle  to  interfere 
with  the  excited  system  which  can  no  longer  be  considered  to 
exist  as  a  free  one.  In  this  respect  our  analysis  is  more  accu¬ 
rate. 

Figure  3  shows  the  dispersion  curves  of  the  free  surface- 
plasma  waves  in  the  absence  of  the  prism  loading.  The 
dashed  and  continuous  curves  show,  respectively,  the  imagi¬ 
nary  part  of  3  and  the  metal  thickness  £  as  a  function  of  the 
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Fig.  4(a)  Reflectivity  (continuous  line)  aa  a  (unction  of  the  incident 
angle  and  the  ratio  of  the  abeolute  square  of  the  magnetic  fields  at  the 
metal  surface  and  inside  the  prism  (dashed  line)  for  t  -  70  nm.  (b) 
Same  as  (a)  but  with  f  ■  10  nm.  Note  the  enhancement  of  the  optical 
fields  relative  to  the  thick  metal  case. 


real  part  of  0.  One  observes  the  divergence  of  t  for  a  value  of 
0  which  corresponds  to  the  degenerate  case  discussed  before. 
For  smaller  values  of  t  one  obtains  two  values  for  both  the  real 
and  imaginary  parts  of  the  propagation  constant,  each  pair 
of  real  and  imaginary  values  corresponding  to  either  the 
short-range  or  the  long-range  mode.  In  the  presence  of  the 
prism  loading  (with  h  »  0.1)  the  Lorentzian  width  shown  in 
Fig.  2  is  approximately  a  factor  of  2  larger  than  the  one  com¬ 
puted  from  lm(0)  in  Fig.  3.  The  chosen  value  of  h,  therefore, 
has  the  effect  of  roughly  equalizing  the  plasma-wave  losses 
associated  with  reradiation  and  dissipation  in  the  metal. 

In  Figs.  4(a)  and  (b)  we  show  the  field  enhancement  in  the 
vicinity  of  the  metal  surface  for  the  parameters  listed  above. 
In  particular,  in  each  figure  the  continuous  line  is  the  ampli¬ 
tude  reflectivity,  and  the  dashed  line  is  the  ratio  of  the  abso¬ 
lute  square  of  the  magnetic  fields  at  the  metal  surface  and 
inside  the  prism  as  a  function  of  the  angle  of  incidence.  Figure 
4(a)  corresponds  to  £  *  70  nm  and  Fig.  4(b)  to  t  =  10  nm  with 
only  the  long-range  mode  shown  in  the  latter  figure.  For  t  * 
10  nm  the  field  enhancement  can  be  seen  to  be  ~  l  order  of 
magnitude  larger  than  that  obtainable  in  the  case  of  the 
thicker  metal  film,  which  corresponds  to  the  Otto  geometry. 
Our  calculations  with  the  Kretschmann  geometry  show  the 
enhancement  to  be  essentially  the  same  as  those  obtained  in 
the  Otto  geometry. 

In  summary,  we  have  presented  a  comparison  between  the 
properties  of  the  surface-plasma  waves  excited  by  three  dif¬ 
ferent  prism-coupling  geometries.  The  analysis  shows  that 
our  geometry  can  give  rise  to  a  field  enhancement  that  can  be 
1  order  of  magnitude  larger  than  those  obtainable  by  the  two 
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the  light  is  guided  along  a  very  thin  metal  film,  for  nonlinear 
interactions  such  as  second  harmonic  generation,  four-wave 
mixing,  and  other  power-dependent  refrcctive-index  phe¬ 
nomena.® 

This  work  is  supported  in  part  by  NSF  grant  ECS-81. 20348. 
One  of  us  (RTD)  -wishes  to  thank  Hyatt  Gibbe  for  an  invitation 
to  spend  a  sabbatical  at  the  Optical  Sciences  Center  he  is  on 
sabbatical  leave  from  the  University  of  Toledo. 


•  ««*•«•«  «MW«« 

1.  A.  Otto,  Z.  Phys.  216, 398  (1968). 

2.  E.  Kretschmann,  Z.  Phys.  241,313  (1971). 

3.  W.  H.  Weber  md  G.  W.  Ford,  Opt.  Lett.  6. 122  (1981). 

4.  D.  Sand.  Pbys.  Rev.  Lett.  47, 1927  (1981). 

5.  M.  Bom  end  E.  Wolf,  Principle  of  Optict  (Pergamon,  New  York, 
1964). 

6.  D.  Sarid,  R.  T.  Deck,  and  J.  J.  Fasano  (submitted  to  J.  Opt.  Soc. 
Am.). 


OUT  IN 


Optical  field  •nhancamwnt  by  long-range  surface- 
plasma  waves 

Dror  Sand,  Robert  T.  Deck,  Alan  E.  Craig,  Robert  K. 
Hickemeil,  Ralph  S.  Jameson,  and  Joseph  J.  Fasano 
Joseph  Fasano  ia  with  U.S.  Amy,  ARRADCOM,  Dover, 
New  Jersey  07801;  the  other  authors  are  with  University 
of  Arizona,  Optical  Sciences  Center,  Tucson,  Arizona 
95721. 

Received  21  May  1982. 

Consider  a  TM  polarized  monochromatic  wave  that  ia  in¬ 
cident  on  the  base  of  a  prism  from  which  it  ia  totally  internally 
reflected.  If  a  metallic  substrate  is  placed  in  the  vicinity  of 
the  prism  so  that  the  evanescent  fields  of  the  reflected  light 
penetrate  the  metal,  and  the  incident  angle  is  appropriate,  a 
surface-plasma  wave  can  be  excited  along  the  metal  gap. 
Otto1  has  analyzed  this  geometry  and  has  shown  both  theo¬ 
retically  and  experimentally  that  by  measuring  the  reflectivity 
of  a  plane  wave  as  a  function  of  the  angle  of  incidence  one 
obtains  both  the  speed  and  the  range  of  the  surface-plasma 
wave  from  the  position  of  the  minimum  and  the  width  of  the 
Lorentzian  absorption  curve,  respectively.  In  his  geometry 
the  metal  thickness  was  much  larger  than  the  decay  length  of 
the  surface- plasma  field  in  the  metaL  Kretschmann2  con¬ 
sidered  a  different  geometry  in  which  a  thin  metal  film  ia  de¬ 
posited  on  the  bottom  of  the  prism  and  a  piano  wave  incident 
on  the  base  of  the  prism  excites  a  surface- plasma  wave  on  the 
lower  side  of  the  metal  film.  The  excitation  in  this  case  takes 
place  via  the  optical  fields  that  penetrate  the  metaL  As  with 
the  Otto  geometry,  the  excitation  is  possible  because  the  speed 
of  the  surface-plasma  wave  is  slower  than  the  speed  of  the  light 
in  the  prism. 

Weber  and  Ford3  have  recently  shown  that  the  fields  ac¬ 
companying  the  surface- plasma  wave  are  larger  than  those 
associated  with  total  internal  reflection  in  the  absence  of 
surface- plasma  waves.  In  their  theory  they  assume  a  general 
coupling  mechanism  that  is  100%  efficient  and  equate  the 
optical  power  density  injected  by  the  coupler  under  steady 
state  conditions  to  the  power  dissipated  by  the  surface-plasma 
wave. 

They  show  that  the  excitation  of  a  surface- plasma  wave  at 
an  optical  frequency  can  give  rise  to  an  enhancement  of  the 
optical  electric  fields  at  a  silver  metal  surface  by  a  factor  as 
large  as  60.  Their  theory,  however,  is  approximate  because 
in  calculating  the  enhanced  fields  they  use  the  free  surface- 
plasma  fields  instead  of  the  true  fields  in  the  presence  of  the 
excitation  medium  (a  prism  or  a  grating,  for  example.) 

In  this  Letter  we  wish  to  show  theoretically  that  (1)  it  is 
possible  to  obtain  an  exact  solution  to  the  field  enhancement 
for  both  the  Otto  and  the  Kretschmann  geometries  and  (2)  our 
new  geometry4  makes  it  possible  to  obtain  field  enhancements 
that  are  as  much  as  an  order  of  magnitude  larger  than  those 
obtainable  with  the  two  other  geometries. 

The  three  prism-coupling  geometries  are  shown  in  Fig.  1, 
where  I  and  III  correspond  to  the  Otto  and  Kretschmann 
geometries,  respectively,  and  II  to  our  geometry.  In  each  case 
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Fig.  1.  Three  prism -coupling  geometries  for  era  ting  e  eurfece- 
plagma  wave:  t,  ID,  end  II  cone* pood  to  the  Otto,  Kretschmann,  and 
our  geometries,  reepectively. 

t  denotes  the  metal  thickness,  and  the  black  dot  in  each  di¬ 
agram  indicates  the  point  where  the  optical  fields  are  evalu¬ 
ated.  In  geometries  I  and  II  the  evanescent  fields  created  at 
the  prism  surface  need  to  penetrate  an  intermediate  medium 
separating  the  prism  and  the  metal  film.  The  exponential 
decay  of  the  field  amplitude  in  this  medium  is  expressed  by 
the  factor  h  *  exp(-KCS),  where  K  m  2x/X,  X  ia  the  optical 
wavelength  in  free  space,  S  and  n?  are  the  thickness  and  re¬ 
fractive  index  of  the  gap,  Cr  *  /J2  -  and  K0  is  the  propa¬ 
gation  constant  of  the  surface-plasma  wave  which  has  only  a 
real  part  for  the  plane  wave  approximation.  In  this  work  we 
assume  a  varying  thickness  5  so  that  h  ■  0.1  for  any  given 
angle.  The  reason  for  this  choice  is  that  under  these  condi¬ 
tions  the  reflectivity  obtains  its  minimum  value  for  the  ap¬ 
propriate  angle  for  most  of  the  cases  we  consider.  It  can  be 
seen  that  geometry  II  merges  into  the  Otto  geometry  for  t 
sufficiently  large  and  h  <  1  and  into  the  Kretschmann  ge¬ 
ometry  for  t  sufficiently  small  and  h  *  1.  For  geometry  II  to 
b«  maximally  effective,  the  refractive  indices  of  the  media 
bounding  the  metal  film  on  either  side  need  to  be  approxi¬ 
mately  the  same. 

Now  in  the  general  case  of  a  stratified  medium  characterized 
by  complex  dielectric  constants  the  reflectivity  of  a  TM 
polarized  monochromatic  plane  wave  at  the  surface  of  the 
medium  and  the  fields  excited  in  each  of  the  media  can  be 
computed  by  standard  techniques.5  For  the  plane  wave  as¬ 
sumption  to  be  applicable  to  our  geometry  it  is  necessary  that 
the  range  of  the  surface-plasma  wave  be  appreciably  smaller 
than  the  width  of  the  incident  optical  beam.  In  the  visible 
range  of  the  spectrum,  where  the  latter  range  is  <1  mm,  this 
requirement  is  easily  satisfied,  and  the  simple  plane  wave 
calculations  yield  accurate  results  for  all  three  geometries. 

To  demonstrate  the  field  enhancement  obtainable  with  our 
geometry  we  choose  parameters  appropriate  to  the  case  of  a 
silver  metal  film  having  a  dielectric  constant  given  by  -  L6  + 
(0.52  at  a  wavelength  of  632.8  nm,  bounded  by  two  dielectric 
media  having  a  refractive  index  of  1.5.  Our  prism  is  assumed 
to  be  rutile  with  a  refractive  index  of  2.875.  Figure  2  shows 
the  amplitude  reflectivity  at  the  base  of  the  prism  for  three 
metal  film  thicknesses:  t  *  2000, 70,  and  50  nm  with  h  *0.1. 
Note  how  the  single  Lorentzian  absorption  curve  associated 
with  the  thick  metal  film  splits  into  two  Lorentzians  as  the 
film  thickness  decreases.  The  splitting  is  a  direct  result  of  the 
lifting  of  the  degeneracy  in  the  propagation  constants  of  the 
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Fig.  2.  FUflvcUvity  «a  a  function  of  tba  incident  angia  inside  tha 
priam  for  our  geometry  where  no  “  nj.  ( — ),  (•••),  and  (— -)  are 
for  t  «  2000, 70,  and  50  no,  respectively.  Nota  tha  splitting  of  tha 
single  Lorentzian  into  two  Lorentzian*  as  tha  metal  film  thickness 
decreases.  The  wide  and  the  narrow  Lorentnans  correspond  to  tha 
short-range  and  long-range  surface-plasma  waves,  respectively. 


Rm (b mba) 

Fig.  3.  Dispersion  curves  of  the  free  surface- plasma  waves:  lm(0) 
(dashed  line)  and  t  (continuous  line)  as  a  function  of  Re<£).  Nota  tha 
divergence  oil  for  Reiff)  that  corresponds  to  tha  thick  metal  case,  and 
trie  two  values  of  Re</3)  and  Im(/3)  for  each  t  for  tha  thinner  metal 

case. 


two  independent  surface-plasma  waves  propagating  on  each 
surface  of  the  thick  metal.  Once  the  metal  film  thickness 
becomes  small  enough,  the  fields  of  these  two  surface- plasma 
waves  inside  the  metal  overlap,  and  the  degeneracy  is  lifted. 
For  such  a  thickness,  one  long-range  and  one  short-range 
surface-plasma  wave  can  be  excited  as  a  single  eigenmode  that 
satisfies  all  the  boundary  conditions  simultaneously.4  One 
observes  that  the  Lorentzian  width  at  half  maximum  for  the 
3hort-range  and  for  the  long-range  mode  increases  and  de¬ 
creases,  respectively.  The  width  consists  of  two  contributions, 
one  from  the  power  dissipation  in  the  metal  film  and  one  from 
the  priam  loading.  This  analysis  differs  from  that  of  Weber 
and  Ford,  because  i .  takes  into  account  the  important  fact  that 
any  efficient  means  of  excitation  has  in  principle  to  interfere 
with  the  excited  system  which  can  no  longer  be  considered  to 
exist  as  a  free  one.  In  this  respect  our  analysis  is  more  accu¬ 
rate. 

Figure  3  shows  the  dispersion  curves  of  the  free  surface- 
plasma  waves  in  the  absence  of  the  prism  loading.  The 
dashed  and  continuous  curves  show,  respectively,  the  imagi¬ 
nary  part  of  3  and  the  metal  thickness  t  as  a  function  of  the 


Fig.  4(a)  Reflectivity  (continuous  line)  ea  a  function  of  tha  incident 
angle  and  tba  ratio  of  tha  absolute  square  of  tha  magnetic  fields  at  tha 
metal  surface  and  inaida  tha  priam  (dashed  line)  for  t  -  70  am.  (b) 
Same  sa  (a)  but  with  t  ■  10  nm.  Nota  the  enhancement  of  tha  optical 
fields  relative  to  the  thick  metal  case. 

real  part  of  0.  One  observes  the  divergence  of  t  for  a  value  of 
8  which  corresponds  to  the  degenerate  case  discussed  before. 
For  smaller  values  of  £  one  obtains  two  values  for  both  the  real 
and  imaginary  parts  of  the  propagation  constant,  each  pair 
of  real  and  imaginary  values  corresponding  to  either  the 
short-range  or  the  long-range  mode.  In  the  presence  of  the 
prism  loading  (with  h  «  0.1)  the  Lorentzian  width  shown  in 
Fig.  2  is  approximately  a  factor  of  2 larger  than  the  one  com¬ 
puted  from  Im(/3)  in  Fig.  3.  The  chosen  value  of  h,  therefore, 
has  the  effect  of  roughly  equalizing  the  plasma-wave  losses 
associated  with  reradiation  and  dissipation  in  the  metaL 

In  Figs.  4(a)  and  (b)  we  show  the  field  enhancement  in  the 
vicinity  of  the  metal  surface  for  the  parameters  listed  above. 
In  particular,  in  each  figure  the  continuous  line  is  the  ampli¬ 
tude  reflectivity,  and  the  dashed  line  is  the  ratio  of  the  abso¬ 
lute  square  of  the  magnetic  fields  at  the  metal  surface  and 
inside  the  prism  as  a  function  of  the  angle  of  incidence.  Figure 
4(a)  corresponds  to  t  »  70  nm  and  Fig.  4(b)  to  t  »  10  nm  with 
only  the  long-range  mode  shown  in  the  latter  figure.  For  t  « 
10  nm  the  field  enhancement  can  be  seen  to  be  M  order  of 
magnitude  larger  than  that  obtainable  in  the  case  of  the 
thicker  metal  film,  which  corresponds  to  the  Otto  geometry. 
Our  calculations  with  the  Kretschmann  geometry  show  the 
enhancement  to  be  essentially  the  same  as  those  obtained  in 
the  Otto  geometry. 

In  summary,  we  have  presented  a  comparison  between  the 
properties  of  the  surface-plasma  waves  excited  by  three  dif¬ 
ferent  prism-coupling  geometries.  The  analysis  shows  that 
our  geometry  can  give  rise  to  a  field  enhancement  that  can  be 
l  order  of  magnitude  larger  than  those  obtainable  by  the  two 
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interactions  such  as  second  harmonic  generation,  four-wave 
mixing,  and  other  power-dependent  refractive-index  phe¬ 
nomena.* 
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The  propagation  and  attenuation  constants  of  long-range  surface- platmon  polariton*  propagating  on  thin  metal 
film*  have  been  measured  for  several  film  thicknesses  by  using  ■  hemispherical  retroreflecting  coupler.  The  results 
are  in  agreement  with  our  theoretical  dispersion  curves,  which  predict  that  the  attenuation  constant  of  this  mode 
can  decrease  by  a  factor  of  63  relative  to  that  obtained  with  a  thick  metal  film. 


It  is  well  known  by  now  that  surface-plasmon  polari-  lens  preceding  the  coupler  ensures  that  the  beam  inside 

tons  that  propagate  on  a  thin  metal  film  split  into  two  it  is  collimated.  The  beam  is  retroreflected  by  the 

branches,  denoted  as  the  short-  and  the  long-range  hemispherical  coupler,  first  from  the  metallized  surface, 

modes,  when  the  thickness  of  the  metal  film  on  which  then  from  the  coupling  surface,  so  that  at  any  given 

they  propagate  decreases  below  a  critical  value  of  about  angle  of  incidence  the  reflected  beam  is  parallel  to  the 

100  nm. 1,2  The  main  features  of  the  long-range  mode,  incident  one.  These  two  beams  are  then  separated  by 

which  is  the  one  of  interest,  are  that,  as  the  thickness  of  a  pellicle  beam  splitter,  and  the  reflected  one  is  focused 

the  metal  film  decreases,  both  its  propagation  and  decay  onto  a  photodetector. 

constants  (the  inverse  of  the  range)  decrease,  while  The  hemispherical  coupler  presses  on  a  thin  metal 

amplitudes  of  the  electric  and  magnetic  fields  associated  film  deposited  on  top  of  a  specially  polished  flat  glass 

with  this  mode  increase.  As  a  result,  nonlinear  inter-  slide.  An  index-matching  fluid  fills  the  gap  separating 

actions  involving  the  long-range  mode,  such  as  sec-  the  coupler  and  the  metal  film,  providing  a  symmetric 

ond-harmonic  generation3  and  the  intensity-dependent  refractive-index  structure  bounding  the  metal  film, 

propagation  constant,4-5  will  be  greatly  enhanced.  The  glass  slide  is  supported  by  three  ball  bearings, 

Previous  measurements  of  the  real  and  imaginary  which  are  spring  mounted  on  separate  micrometers, 

components  of  the  propagation  constant  of  surface-  The  gap  can  therefore  be  adjusted  in  both  magnitude 

plasmon  polaritons  propagating  on  thin  metal  films  and  shape  to  conform  to  the  optimal  coupling  specifi- 

(down  to  50  nm)  have  indicated  that  both  of  these  cations. 

components  decrease  as  the  metal  film  thickness  de-  The  microcomputer  drives  the  stepping  motor,  which 

creases.6  In  this  Letter  we  present  experimental  results 
in  which  the  reflectivity  from  very  thin  metal  films 
(down  to  10  nm)  was  measured  as  a  function  of  the  angle 
of  incidence.  These  results,  obtained  by  a  computer¬ 
ized  system  that  incorporates  a  hemispherical  retrore¬ 
flecting  coupler  together  with  an  index-matching  fluid, 
are  in  agreement  with  our  previously  published  theo¬ 
retical  dispersion  curves. 

In  fitting  the  reflectivity  measurements  with  the 
theory,  we  use  only  the  Fresnel  equations  with  complex 
refractive  indices  and  the  dispersion  equation  of  the  free 
surface-plasmon  polaritons. 

Our  system  consists  of  three  parts,  the  computer- 
controlled  electronics,  the  optical-beam-shaping  com¬ 
ponents,  and  the  hemispherical  retroreflecting  coupler. 

We  use  a  microcomputer  to  drive  a  step-motor-driven 
rotary  table  with  a  reproducible  angular  resolution  of 
1/240  deg. 

Light  from  a  He-Ne  laser  operating  at  a  wavelength 
of  632.8  nm  passes  through  a  1-kHz  chopper,  a  spatial  Fig  L  The  hemi8pherical  retroreflecting  coupler  showing 

filter,  a  beam  expander,  a  square  aperture,  and  a  the  three  ball-bearing  pressure  points  for  shaping  the  gap  that 

beam-shaping  lens.  The  beam  is  next  incident  upon  a  separates  the  coupler  from  the  metal  film.  An  index- 

hemispherical  coupler  that  was  cut  in  half,  metallized,  matching  fluid  fills  the  gap  and  makes  it  possible  to  excite  the 

and  then  glued  back  together,  as  shown  in  Fig.  1.  The  long-range  surface-plasmon  polariton. 
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in  turn  rotates  the  turntable,  while  the  detector  moni¬ 
tors  the  intensity  of  the  reflected  beam.  The  signal 
from  the  photodetector  is  fed  into  a  lock-in  amplifier, 
converted  into  a  digital  signal,  and  transferred  to  the 
microcomputer.  A  plotter  then  records  the  reflected 
beam  intensity  as  a  function  of  the  beam  angle  of  inci¬ 
dence  on  the  base  of  the  coupler. 

In  the  vicinity  of  the  resonance  angle  of  the  surface- 
plasmon  polariton,  the  reflected  beam  decreases  in  in¬ 
tensity  because  of  the  dissipation  of  the  optical  power 
carried  by  the  excited  mode.  The  FWHM  of  the  reg¬ 
istered  Lorentzian  absorption,  together  with  the  angle 
of  maximum  absorption,  is  measured.  The  pressure  on 
each  of  the  ball  bearings  is  adjusted  several  times,  and 
for  each  adjustment  the  width  and  the  location  of  the 
Lorentzian  absorption  are  measured.  From  these  re¬ 
sults  one  obtains  the  natural  width  and  location  of  the 
Lorentzian  absorption  that  is  due  to  the  dissipation  of 
the  optical  power  by  the  surface-plasmon  polariton. 
The  real  and  imaginary  components  of  the  dielectric 
constant  are  determined  by  adjusting  their  value  in  the 
theoretical  reflectance  calculation  so  that  the  angular 
location  and  width  of  the  resonance  match  those  of  the 
experimental  curve  for  the  200-nm-thick  metal  film. 
Once  the  dielectric  constant  is  known,  it  is  used  to  plot 
the  real  and  imaginary  components  of  the  propagation 
constant  for  the  two  branches.  The  dielectric  constant 
that  we  measured  for  the  thick  film,  and  that  we  have 
used  for  the  analysis  of  the  thinner  films,  is  -18.9  + 
i'0.581,  which  is  typical  of  Ag. 

For  thinner  films,  we  monitor  their  nominal  thickness 
during  evaporation  and  then  adjust  the  theoretical 
thickness  to  a  value  such  that  the  theoretically  plotted 
Lorentzian  absorption  has  its  minimum  at  the  same 
angle  as  that  of  the  experimental  one.  We  than  take  the 
ratio  of  the  widths  of  the  experimental  and  the  theo¬ 
retical  Lorentzian  absorptions  and  multiply  the  theo¬ 
retical  imaginary  component  of  the  propagation  con¬ 
stant  by  the  same  ratio  to  obtain  the  experimental 
value. 

We  tested  four  Ag  films,  which  were  vacuum  depos¬ 
ited  slowly  enough  to  permit  approximate  thickness 
measurement.  For  films  thinner  than  30  nm,  the 
thickness  above  which  Ag  agglomeration  islands  dis¬ 
appear,  the  substrate  was  electron  bombarded  during 
evaporation.  The  films  had  the  nominal  thicknesses 
of  2000,  50,  13,  and  10  nm.  Each  film  thickness  was 
monitored  to  an  accuracy  of  about  20%  at  evaporation. 
The  refractive  indices  of  the  coupler,  the  oil,  and  the 
glass  were  1.7788, 1.50173,  and  1.50231,  respectively. 

The  solid  and  dotted  curves  of  Fig.  2  show  the  theo¬ 
retical  values  of  the  real  and  the  imaginary  components, 
respectively,  of  the  propagation  constant  of  the  two 
surface-plasmon  polariton  modes.  One  observes  that, 
as  the  metal  film  thickness  decreases,  both  the  real  and 
the  imaginary  components  decrease  for  one  mode  and 
increase  for  the  other  one,  which  makes  it  possible  to 
distinguish  between  the  two  experimentally.  The 
agreement  between  the  theoretical  and  the  experi¬ 
mental  imaginary  components  is  only  within  a  factor  of 
2  because  of  differences  in  the  effective  dielectric  con¬ 
stants  of  the  thick  and  thin  films. 

The  trend  toward  lower  losses  (and  concomitant 
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Fig.  2.  The  theoretical  results  of  the  real  (dotted  curve)  and 
the  imaginary  (solid  curve)  components  of  the  propagation 
constant  as  a  function  of  the  metal  film  thickness  t.  The 
upper  and  lower  branches  describe  the  short-  and  long-range 
surface-plasmon  polaritons,  respectively.  The  experimental 
values  of  the  real  component  of  the  propagation  constant  are 
marked  by  open  ovals.  The  imaginary  component  ( X )  agrees 
within  a  factor  of  2  with  the  theoretical  results. 


longer  propagation  distances)  for  modes  propagating 
on  thinner  metal  films  is  clear.  For  the  case  of  the 
10-nm-thick  silver  film,  the  thinnest  film  for  which 
measurements  of  the  angular  width  of  the  resonance 
have  been  made,  0.055  deg  at  FWHM  was  recorded.  A 
substantial  portion  of  this  angular  width  is  due  to  an¬ 
gular  divergence  of  the  laser  beam  at  the  coupling  in¬ 
terface,  estimated  to  be  0.044  deg.  The  difference  be¬ 
tween  these  two  angles  is  contributed  by  the  surface- 
plasmon  polariton,  from  which  we  find  that  its  range  is 
0.25  mm.  This  value  is  63  times  larger  than  the  corre¬ 
sponding  value  when  a  thick  metal  film  is  used. 

We  have  also  observed  the  short-range  mode.  It 
exhibits  a  Lorentzian  absorption  having  a  FWHM  that 
broadens  and  shifts  farther  from  the  critical  angle  as  the 
thickness  of  the  film  is  progressively  decreased. 

The  procedure  used  to  compare  the  experimental 
results  with  the  theory  neglects  the  important  contri¬ 
bution  that  the  preparation  method  of  the  film,  the 
creation  of  cermets,  surface  roughness,  and  aging  heve 
on  the  value  of  the  dielectric  constant.  Nevertheless, 
the  experimental  results  indicate  clearly  that  the  real 
and  the  imaginary  components  of  the  propagation 
constant  follow  the  theoretical  prediction  and  that  films 
as  thin  as  10  nm  can  be  readily  measured.  Since  the 
resonances  obtained  for  the  thin  films  are  narrow,  ac¬ 
curate  measurements  of  small  imperfections  can  easily 
be  identified. 

In  summary,  we  present  experimental  results  of  the 
propagation  and  decay  constants  of  the  long-range 
surface-plasmon  polariton  and  show  that  they  are  in 
agreement  with  the  theoretical  dispersion  curves.  The 
measurements  were  carried  out  down  to  10-nm  thin 
metal  films,  for  which  the  range  is  increased  by  a  factor 
of  63  relative  to  that  obtained  with  thick  films. 


In  the  process  of  writing  this  Letter  we  have  learned 
of  the  independent  measurement  of  the  long-range 
surface-plasmon  polariton  by  J.  C.  Quail  et  al.  of  the 
University  of  Toledo,  who  were  kind  enough  to  send  us 
a  preprint  of  their  work.  The  authors  would  like  to 
thank  Steve  Browning  and  Fred  Van  Milligen  for  the 
preparation  of  the  films.  This  work  is  supported  by 
National  Science  Foundation  grant  ECS-8120348,  De¬ 
fense  Advanced  Research  Projects  Agency  grant 
N60530-83-C-0089,  and  by  the  U.S.  Air  Force  Office  of 
Scientific  Research  (AFSC)  and  the  U.S.  Army  Re¬ 
search  Office. 
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Coupling  between  finite  electromagnetic  beam  and  long-range 
surface-plasmon  mode 


Robert  T.  Deck,  Dror  Sarid,  Grieg  A.  Olson,  and  J.  M.  Eison 


In  this  paper,  we  consider  a  layered  structure  consisting  of  dielectric,  metal,  and  dielectric  media  and  com¬ 
pute  the  fields  produced  in  the  structure  by  a  finite  width  beam  of  incident  radiation  which  couple*  to  a  long- 
range  surface-plasmon  mode  in  the  metallic  layer.  The  analysis  allows  for  complex  dielectric  constants  in 
all  media  and  focuses  on  optical  radiation.  Curves  are  presented  which  show  the  profiles  of  the  reflected 
field  intensity  for  different  values  of  the  propagation  length  of  the  mode  relative  to  the  width  of  the  incident 
beam.  We  discuss  the  conditions  in  which  the  field  reradiated  from  the  long-range  surface  mode  can  inter¬ 
fere  with  the  specularly  reflected  radiation  so  at  to  produce  an  interference  zero  in  the  reflected  profile.  In 
these  conditions  our  analysis  allows  for  a  determination  of  the  separate  loaa  coefficients  of  the  mod*  associ¬ 
ated  with  energy  reradiation  and  dissipation,  respectively 


I.  Introduction 

There  is  continued  interest  in  application  of  elec¬ 
tromagnetic  energy  propagation  along  boundary  in¬ 
terfaces  and  in  thin  film  waveguides.  The  modes  of 
field  propagation  involved  will  be  referred  to  here  as 
boundary  layer  modes.  The  fields  of  such  modes  ef¬ 
fectively  propagate  on  or  parallel  to  the  interface  sur¬ 
faces  of  a  layered  structure  consisting  of  media  of  dif¬ 
ferent  dielectric  properties. 1  Because  the  fields  of  the 
modes  do  not  couple  directly  to  volume  electromagnetic 
fields,  the  excitation  and  detection  of  the  modes  require 
certain  techniques. 

A  general  description  of  the  modes  in  the  case  of 
lossless  thin  film  dielectric  waveguides  has  been  given 
by  Tien  and  Ulrich-  and  Ulrich.1  The  present  paper 
presents  a  similar  description  in  the  case  where  the 
guiding  film  is  metallic  and  the  media  which  bound  it 
are  in  general  lossy.  In  this  case  the  propagating 
bound -,ry  layer  modes  correspond  to  surface-plasmon 
modes  on  the  metal  film. 

A  large  number  of  papers  have  discussed  the  excita¬ 
tion  and  detection  of  surface-plasmon  modes  on  the 
surface  of  a  thick  metal  film  (e.g.,  Refs.  4  and  5). 
However,  recent  work6  has  pointed  to  a  new  type  of 
surface-plasmon  mode  that  can  propagate  on  the  sur- 
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faces  of  a  thin  metal  film.  This  new  mode  referred  to 
as  a  long-range  surface-plasmon  wave  (LRSPW)  has  a 
range  of  propagation  which  can  be  more  than  an  order 
of  magnitude  larger  than  ranges  previously  observed. 
Because  this  new  mode  has  possible  applications  as  a 
sensitive  probe  of  surface  roughness  and  the  dielectric 
properties  of  materials,  we  feel  that  it  is  important  to 
have  an  overall  description  of  the  fields  associated  with 
such  modes  and  of  the  effects  which  their  presence  has 
on  the  reflectivity  of  ideal  multilayered  surfaces.  We 
give  here  a  description  of  both  in  the  case  where  the 
modes  are  excited  through  a  prism  coupler.  The  re¬ 
flectivity  is  calculated  in  the  cases  where  the  incident 
beam  has  both  an  infinite  and  a  finite  width.  In  the 
latter  case,  with  the  incident  beam  assumed  to  be  a 
square  wave,  we  describe  the  profile  of  the  reflected 
beam  in  the  presence  of  a  long-range  surface-plasmon 
mode. 

The  effects  of  boundary  layer  modes  on  the  profiles 
of  reflected  beams  have  been  described  in  previous 
papers  (e.g.,  Refs.  3,  7,  and  8)  in  which  losses  due  to 
dissipation  are  for  the  most  part  neglected.  A  de¬ 
scription  of  the  reflected  profile  of  a  Gaussian  incident 
beam  in  the  case  of  surface-pk  ?mon  excitation  in 
a  thick  metal  film  is  given  in  Ref.  9.  Because  the 
propagation  length  of  a  surface  plasmon  wave  on  a  thick 
metal  film  decreases  rapidly  as  the  frequency  of  the 
wave  increases,  observable  effects  are  obtained  in  Ref. 
9  only  at  frequencies  in  the  IR.  In  the  present  paper 
where  we  focus  on  the  excitation  of  a  LRSPW  on  a  thin 
metal  film,  observable  effects  can  be  obtained  even  at 
frequencies  corresponding  to  visible  radiation,  and  in 
our  applications  we  consider  such  radiation  exclu¬ 
sively. 
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Fig.  1.  (a)  Geometry  of  boundaries  between  three  media  which  de¬ 
fine  iayered  structure.  <  b»  Geometry  used  to  excite  boundary  layer 
mode  on  layered  structure  2-1-0. 

In  Sec.  II. A  we  present  the  general  theory  of  excita¬ 
tion  of  a  surface-plasmon  mode  by  a  plane  electro¬ 
magnetic  wave,  with  the  focus  on  the  case  where  the 
coupling  between  the  incident  and  plasmon  fields  can 
be  considered  weak.  By  use  of  a  simple  pole  approxi¬ 
mation  in  this  case  we  develop  simplified  expressions 
for  the  reflection  and  transmission  coefficients  in  the 
presence  of  the  plasmon  mode.  Here  the  notation  and 
style  of  Ref.  3  are  followed,  and  for  completeness  certain 
details  given  there  are  repeated.  In  Sec.  II.B  the  plane 
wave  theory  is  made  use  of  to  develop  expressions  for 
the  reflection  and  transmission  coefficients  in  the  case 
of  plasmon  excitation  by  a  finite  width  electromagnetic 
wave.  The  results  we  obtain  represent  generalizations 
of  the  equations  of  Ref.  3  for  the  case  of  lossy  media.  In 
particular,  we  allow  for  complex  dielectric  constants  in 
all  media  except  the  coupling  prism.  We  make  use  of 
our  results  applied  to  the  case  of  an  incident  square 
wave  to  describe  a  method  of  analysis  by  which  it  is 
possible  to  determine  the  separate  loss  coefficients  of 
the  mode  associated  with  energy  reradiation  and  dis¬ 
sipation. 

In  Sec.  Ill  the  formulas  developed  in  Sec.  II  are  ap¬ 
plied  to  a  specific  case  of  excitation  of  a  long  range 
surface-plasmon  mode  by  a  square  wave  He-Ne  beam. 
Here  we  present  curves  showing  the  profiles  of  the  re¬ 
flected  field  intensity  for  various  ranges  of  propagation 
of  the  mode  relative  to  the  width  of  the  incident  beam. 
Finally  in  the  Appendix  we  give  general  formulas  for  two 
quantities  associated  with  the  surface  modes  consid¬ 
ered. 

II.  G«n«ral  Theory 

A.  Mode  Excitation  by  a  Plane  Wave 

We  are  interested  in  the  excitation  of  a  boundary 
layer  mode  of  the  electromagnetic  field  within  the 
layered  structure  defined  by  the  2-1  and  1-0  boundary 
surfaces  shown  in  Fig.  1(a).  By  definition  the  fields  of 
such  a  mode  must  decay  exponentially  away  from  both 
the  *.-1  surface  (in  the  +z  direction)  and  the  1-0  surface 


(in  the  —  z  direction)  so  as  to  propagate  effectively  only 
parallel  to  the  boundary  interfaces.  For  such  a  mode 
to  exist  the  component  of  the  propagation  vector  par¬ 
allel  to  the  interfaces  k  must  satisfy  a  so-called  dis¬ 
persion  equation  expressible  in  the  form--' 

-  '.V'i;  -  »  -mw.  "i  ■  ii. I .  ill 

Here  use  is  made  of  the  notation 


1 21 


where  r,/  is  equal  to  the  amplitude  reflection  coefficient 
for  reflection  of  a  plane  wave  of  wave  vector  k.  *  (k  x 
+  kj:£)  at  a  boundary  surface  separating  two  media  with 
dielectric  constants  t,  and  «/,  respectively.  In  the  case 
(of  interest)  of  a  TM  polarized  wave 


r,i  -  exp(-2n»,i)  « 


The  dispersion  equation  (1)  can  be  derived  from  the 
(unphysical)  requirement  that  the  total  reflection 
coefficientfor  reflection  from  the  2-1-0  (or  0-1-2)  in¬ 
terface  be  equal  to  infinity.  In  particular,  with  the  total 
reflection  coefficient  for  the  2-1-0  interface  written 
as 


r>i  +r,iiMp(2i*l,d|) 

I  + 

the  condition  for  a  boundary  layer  mode  derives  from 
the  requirement  that  the  denominator  of  r2lo  vanish, 
or  equivalently,  that  a  pole  occur  in  the  quantity  [1  — 
rnrw  exp(2i*lldI)]-1. 

It  is  convenient  to  introduce  the  notation 


ffo(0)  ■  ri2r,0exp(2i'A urfp.  (6) 

With  dm  defined  to  be  a  value  of  8  for  which  Ro  equals 
unity,  the  pole  condition  equivalent  to  Eq.  ( 1 )  assumes 
the  form 


ft..liiw)  *  exp|2i(*ud|  -  iDf,  -  •/>iii)||.i,  *  exp|i<J.iid„)|  -  1.  ill- 

Expressed  as  the  condition  on  the  reflection  coefficient, 
t 210  *  « ,  the  mode  condition  defines  a  boundary  layer 
mode  on  the  2-1-0  interface  as  a  mode  which  persists 
in  the  absence  of  an  incident  field  in  medium  2.  On  the 
other  hand,  since  the  dispersion  relation  (1)  is  found  to 
have  only  solutions  with  kj  greater  than  y/t2  (w/c),  the 
k  -vector  of  the  boundary  layer  mode  must  exceed  the 
propagation  vectors  of  bulk  electromagnetic  waves  in 
medium  2,  and  the  mode  can,  therefore,  be  excited  only 
by  an  evanescent  wave  in  that  medium. 

In  a  so-called  Otto  configuration10  an  evanescent 
wave  with  k  =  (w/c)  \'e3  sin#,-  >  (w/c)  y  <2  is  excited  in 
medium  2  by  an  electromagnetic  wave  which  is  totally 
internally  reflected  at  a  boundary  surface  3-2  that  is 
separated  from  the  2-1  surface  by  a  gap  of  distance  d  > 
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[as  in  Fig.  l(b)j.  The  total  internal  reflection  condition 
requires  <3  >  to.  In  order  that  the  amplitude  of  the 
evanescent  field  be  large  enough  at  the  2-1  surface  to 
excite  an  observable  mode  on  the  2-1-0  interface,  the 
amplitude  factor 


tion  r0i23  =  ®,  which  amounts  to  the  requirement  that 
a  pole  occur  in  the  quantity 

u(0)*|l  -  r)1.;i7-ioexp(2i>i..d]l|"1  (101 

Because  r123  is  equal  to  the  expression 


h  =  expUA^ds)  »  exp 


(T) 


r u  +  exp(2jk  >jdj) 

1  + 


ill) 


must  be  sufficiently  larger  than  zero,  which  requires 
that  do  be  sufficiently  small. 

In  the  case  of  interest  in  which  medium  1  is  metallic, 
the  dielectric  constant  fi  of  medium  1  is  a  complex 
number  with  a  negative  real  part,  and  the  fields  of  a 
boundary  layer  mode  on  the  2-1-0  interface  are  evan¬ 
escent  within  medium  1  as  well  as  within  media  2  and 
0.  In  this  case,  for  TM  polarized  radiation,  there  exist 
for  the  mode  number  m  =  1  in  Eq.  (1)  two  complex  so¬ 
lutions  of  the  dispersion  equation  for  i3m.  If  the 
thickness  d  1  of  the  metallic  medium  is  sufficiently  large, 
these  two  solutions  correspond  to  independent  sur- 
face-plasmon  modes  on  the  separate  surfaces  2-1  and 
1-0.  In  particular,  where 


exp(2j>i,d()  -  exp 


dlv/pPrreT 


is  approximately  zero  in  Eq.  (6)  (for  large  di),  the  dis¬ 
persion  equation  Ro  (/3m)  =  1  reduces  to  the  separate 
equations 


-1.0, 

rl0 

(8a) 

—  »0, 
r\2 

(1*21  m 

(8b) 

which  correspond  to  dispersion  relations  for  indepen¬ 
dent  surface-plasmon  modes  on  the  interfaces  1-0  and 
1-2,  respectively.  As  the  thickness  d  1  of  the  metallic 
medium  is  reduced,  however,  the  evanescent  tails  of 
separate  modes  on  the  two  interfaces  overlap,  and  the 
formerly  independent  modes  interact  to  produce  two 
distinct  modes  characterized  by  electric  field  distribu¬ 
tions  in  the  metal  that  are,  respectively,  symmetric  and 
antisymmetric  with  respect  to  the  transverse  direction. 
Because  the  mode  with  the  antisymmetric  field  distri¬ 
bution  has  a  smaller  fraction  of  its  total  field  inside  the 
metallic  medium  and,  therefore,  suffers  less  dissipation, 
it  has  a  longer  propagation  range  than  does  the  other 
mode.  It  is  this  long-range  mode  which  is  of  primary 
interest  in  the  present  paper.  Since  the  range  of  the 
mode  is  determined  by  the  imaginary  part  of  the  mode 
k -vector,  the  long-range  mode  corresponds  to  a  small 
value  of  lmd. 

In  the  presence  of  the  3-2  boundary  the  reflection 
coefficients  r210  and  r0I2  at  the  upper  and  lower  surfaces 
of  the  layered  structure  need  to  be  replaced  by  the  re¬ 
spective  coefficients  r3 210  and  r0l23  with  r,//m  defined  by 
the  equation 


_  r„  +  r;im  exp(2ikltd,) 

rt/(m  = 

1  +  r„r,im  exp(Zik„d,l 

The  condition  on  d  required  for  a  boundary  layer  mode 
on  the  2-1-0  interface  is  then  expressible  as  the  condi- 


provided  h2  -  exp^ik^do)  is  sufficiently  small,  r  123  » 
r l2  +  0 (h2),  and  the  poles  of  the  quantity  u(fi)  are 
shifted  only  slightly  from  the  poles  of  the  quantity  [1  - 
/?o(/3))— 1  -  [1  -  ri2r10exp(2i7?izd1)]_l.  It  follows  that 
the  values  of  d  corresponding  to  modes  of  the  interface 
2-1-0  in  the  presence  of  the  3-2  boundary  are  ap¬ 
proximately  the  same  as  those  of  the  free  interface. 
The  shifts  in  the  positions  of  the  poles  in  the  presence 
of  the  additional  boundary  result  from  leakage  of  energy 
out  of  the  modes  in  the  form  of  radiation  into  medium 
3.  This  leakage  adds  an  additional  damping  term  to  the 
mode  propagation  vector  which  appears  as  an  additional 
positive  component  of  lmdi- 

1.  Case  of  Weak  Coupling 

Where  the  parameter  h  of  Eq.  (7)  is  large,  so  that 
there  is  a  strong  coupling  between  medium  3  (the  prism) 
and  medium  1  (the  metal),  the  additional  damping  term 
in  the  mode  fields  resulting  from  radiation  into  medium 
3  causes  boundary  layer  modes  on  the  2-1-0  interface 
to  be  rapidly  extinguished.  In  the  cases  of  most  phys¬ 
ical  interest,  therefore,  the  parameter  h  is  restricted 
between  limits  imposed  by  the  dual  requirements  that 
h2  be  sufficiently  greater  than  zero  to  allow  for  excita¬ 
tion  of  an  interface  mode  and  sufficiently  less  than  one 
to  prevent  the  rapid  extinction  of  the  mode.  In  practice 
here  we  restrict  the  values  of  h 2  to  lie  in  the  interval 
between  0.01  and  0.1.  In  this  case  the  coupling  between 
media  1  and  3  is  referred  to  as  weak. 

To  compute  the  reflectivity  of  the  3-2-1-0  layered 
structure  in  the  conditions  of  excitation  of  a  surface- 
plasmon  mode  on  the  2-1-0  interface  we  need  to  eval¬ 
uate  the  amplitude  reflection  coefficient  r32io  for  re¬ 
flection  of  a  plane  wave  field  from  the  3-2  boundary  in 
the  presence  of  the  boundary  layer  mode.  For  this 
purpose  it  is  useful  to  reexpress  r32io  in  terms  of  the  pole 
function  u(&)  of  Eq.  (10)  in  the  form 

^3210  ”  v(j8)u(/3),  (12) 

where 

r3 2  +  ',33r;]rioexp(2i6u<i])  +  ri2[r21  +  rw  exp(2ikud,)j 
(1  +  h2r\2ri3) 

(13) 

Under  the  assumption  of  weak  coupling  with  h  2  «  1, 
the  quantity  r123  in  u($)  can  be  expressed  to  first  order 
in  h2  in  the  form 

rl23  "  r,2(l  +  2ih -r23 sin2«i2).  1 1 4 1 

and  uil3)  can  be  approximated  by  the  expression 

ul/i)  -  (i  -  R(J)|-1.  1101 

where 
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Rid)  »  rl2rwexpl2iklsdiHl  +  2ih  V23  sin2<£i2) 

5K<3>  1 


!  Rotji 


1  +■ 


/?o(3) 

5^(J>  *  2h:r: ,  sin2c/>|  j. 


1  expjiyto(tf)  +  i5^(3)|  >  (15) 


The  poles  of  u(3)  occur  at  values  of  3  denoted  3m  for 
which  #(3)  =  1.  Equivalently  the  positions  3m  of  the 
poles  are  defined  by  the  equation 


l£o(ftn)  +  ”  2mir, 

m-  0,1,2 .  (16) 

with  the  value  of  3m  shifted  from  the  values  3m  corre¬ 
sponding  to  the  modes  of  the  free  2-1-0  interface  by  an 
amount  of  order  h-: 


3m  ®  dm  4*  ^dm- 


(17) 


By  Taylor  expansion  of  and  b\p  about  3m  =  3m  and 
use  of  Eq.  (1),  Eq.  (16)  can  be  reexpressed  as  an  ap¬ 
proximate  equation  for  b$m: 


^"•dm) 


»dm  + 


5^(dm)  2m*’, 


6dm  =  — Ol/'fdm) 

The  imaginary  part  of  53m  determines  the  additional 
damping  of  the  mth  boundary  layer  mode  caused  by  the 
presence  of  the  3-2  boundary.  A  general  expression  is 
given  in  the  Appendix  for  the  quantity 

For  a  value  of  3  in  the  vicinity  of  the  value  3m  the 
function  u(0)  has  the  form  of  a  simple  pole  term  with 
a  coefficient  determined  by  its  Taylor  expansion  about 
3m  The  result  can  be  expressed  in  the  approximate 
form 


2h2r-t3  sin2d>t2 


0m 


(18) 


dm 


u(f})  » 


dd  0- 


(d-3m) 


(19) 


To  the  same  approximation  the  factor  u(3)  of  Eq.  (13) 
can  be  written  for  3  in  the  vicinity  of  3m  as 


<-M(d) 


/3-dm 


5dm 
^^(dm ). 


(20) 


and  the  reflection  coefficient  r32io  can  be  reduced  to  the 
form 


'■3210(d)  ^  r:,2(dm) 


d-dm- 


5dm 


r^(dm)  I 


(d-dm) 


(21) 


where  the  slowly  varying  quantity  r32(3)  in  (20)  and  (21) 
is  approximated  by  its  value  at  3  =  3m- 

Equation  (21)  represents  a  valid  approximation  to 
r  )2io  for  a  range  of  values  of  3  centered  on  3m  and  ex¬ 
tending  on  either  side  of  this  value  about  half  of  the 
distance  to  the  next  pole  of  u  (3).  Since  the  interval 
A3m  between  the  poles  of  u<3)  (in  the  vicinity  of  3  * 
dm )  can  be  computed  from  the  equation 


Allddm)*^!  Adm  “2t.  (22) 

ap\tm 

it  results  that  Eq.  (21)  is  valid  for  a  range  of  3  values 
given  approximately  by 

•^""2r/dd| i.‘  (22)' 

By  use  of  Eq.  (17),  Eq.  (21)  can  be  reexpressed  in  an 
alternate  form  as  the  equation 


r3210(d)  “  ',32(dm) 


1  +■ 


5  (3m 

!_  i 

1 

&  -  +  6&m)  +  ilm(0m  +  60* )J 

(23) 

For  3  in  the  vicinity  of  3m  the  absolute  square  of  the 
bracket  on  the  right-hand  side  of  this  equation  exhibits 
a  Lorentzian  dip  as  a  function  of  3  centered  at  3  = 
Re(3m  +  5 3m).  with  a  width  T  given  by  Im(3m  +  53m) 
=  Im3m-  Therefore,  since  3  is  determined  by  the  angle 
of  incidence  0,  on  the  3-2  boundary  by  the  relation 


d  “  —  ■  —  k  sin#*  »  y/7j  sinBi,  (24) 

(4)  U) 


the  reflectivity  R  defined  by  |r3210(3)|2  exhibits  a  Lo¬ 
rentzian  shape  as  a  function  of  sind,  with  its  center  at 
(0.)o  given  by  y/T^  sin(0,)o  =  Re(3m  +  53m)  =  Re3m- 
Where  the  metallic  medium  1  is  sufficiently  thin  in  this 
case,  the  two  solutions  of  the  mode  dispersion  relation 
form  =  1,  corresponding  to  the  distinct  surface-plasmon 
modes  on  the  metal  film,  appear  in  the  reflectivity  as 
separate  Lorentzians  whose  angular  positions  and 
widths,  respectively,  determine  the  real  and  imaginary 
parts  of  the  3m  values  for  the  separate  modes.  The 
Lorentzian  corresponding  to  the  long-range  surface- 
plasmon  mode  is  distinguished  by  its  narrower  width 
from  that  corresponding  to  the  shorter  range  mode. 

By  way  of  the  approximations  used  to  obtain  the 
above  result  for  r32io  it  is  possible  to  derive  a  reduced 
expression  for  the  amplitude  transmission  coefficient 
1 321  connecting  the  amplitude  of  an  incident  field  in 
medium  3  and  the  amplitude  of  the  resulting  field  in 
medium  1.  In  general  the  resulting  field  in  medium  1 
is  the  resultant  of  the  field  transmitted  through  the  2-1 
boundary  and  the  fields  multiply  reflected  at  the  1-0 
and  1-2-3  interfaces.  Since  the  effect  of  the  latter  re¬ 
flected  fields  on  the  amplitude  of  the  field  in  medium 
l  is  expressible  in  terms  of  the  z  -dependent  factor 

a(Qj)  »  |exp(i/tuz)  +  rl0  exp[-ikx,z  +  2ikXi(dx  +  d2)]i,  (25) 

it  follows  that  the  ratio  of  the  amplitudes  of  the  fields 
in  media  1  and  3  is  given  by  the  product 


a((3^U32i(|(3). 

where,  for  the  geometry  of  Fig.  1,  £32i  can  be  written  in 
terms  of  u(3)  in  the  form 


(.12l($) 


(1  +  '•32)0  +  rai) 

1  +  r 32''2i  exp<2 )*2jd2) 


ulfS)  explil*^  -  *w)d2)  (26) 


Use  of  Eqs.  (19)  and  (25)  and  mode  condition  (1)'  sub¬ 
sequently  reduces  the  product  a  (3, z  >£321  (3)  to  the  ap¬ 
proximate  result 
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B.  Mode  Excitation  by  Finite  Width  Beam 

The  (plane  wave)  results  of  the  previous  section  can 
be  applied  to  the  interaction  of  a  finite  width  electro¬ 
magnetic  beam  with  the  multilayered  structure  of  Fig. 
1(b)  by  decomposition  of  the  incident,  reflected,  and 
refracted  beams  into  plane  wave  components.  Consider 
in  particular  a  beam  with  finite  width  in  the  x-z  plane 
incident  on  the  boundary  surface  3-2  at  z  *  0.  For 
simplicity  we  can  assume  that  the  beam  has  an  unlim¬ 
ited  width  in  the  y  direction  and  consists  of  radiation 
of  a  single  frequency  w.  Since  the  conditions  for  exci¬ 
tation  of  a  surface-plasmon  mode  require  TM  radiation, 
we  take  the  beam  to  be  TM  polarized  with  its  magnetic 
field  vector  perpendicular  to  the  x-z  plane.  Because 
the  electric  field  components  of  the  beam  are  related  to 
the  magnetic  field  component  by  Maxwell’s  equations, 
it  is  sufficient  to  consider  only  the  latter  field  compo¬ 
nent.  At  z  =  0  this  component  has  the  general  form 


By  introduction  of  Fourier  transforms  of  r;tiIo(d)  and 
a(dA)t32i($)  and  use  of  the  convolution  theorem,  Eqs. 
(32)  and  (33)  can  be  rewritten  as  the  equations 

uj/c  c  * 

fl3lx)  «• —  (  dx'RiHuix  -  x')fli(x').  i:)4a) 

2x  J -- 

ui/c  c  * 

Bi(x,z)  ■* —  |  dxTm(x  -  .w)fi3(x'),  (35a) 

2t  J -- 

where 

Rmoix)  *  J*  dfrm o<3)  exp  |i  j  3xJ  •  (34b) 

Tm(xj)  ■  J*  d0a(Qj)t32\(d)  exp  - Bxj  •  (35bt 

Under  the  assumption  that  the  incident  field  B3(x ) 
is  well  collimated,  the  range  Ad  of  d  values  over  which 
(33(d)  is  nonzero  is  determined  strictly  by  the  spatial 
width  w  of  the  incident  beam  through  the  relation 
iu/c)A@w  =*  2ir  or 

A 3  ar  2x/(<i >/c)w.  (36) 

Provided  A  0  is  less  than  or  equal  to  the  approximate 
spacing  (in  d  space)  between  the  surface  modes  of  the 
2-1-0  boundary  interface,  as  defined  by  Eq.  (22)',  the 
function  63(d)  will  overlap  strongly  with  only  one  of  the 
resonance  poles  of  the  functions  r:tj|i,(d)  and  f  121(d)  in 
the  integrals  (32)  and  (33),  and  in  general  the  incident 
beam  will  excite  only  one  of  the  modes  of  the  interface. 
For  (c7w)&'’  in  the  vicinity  of  d,„  the  requirement  Ad  ^ 
Ad,„  imposes  on  «’  only  the  condition 


B3(x.t>  ■  B3y(x^,f)|,-o3  B3(x)  expl-(uf),  (29) 
B3(x)  *  F(x)  exp <ife?x), 

where  6?  represents  the  average  value  of  the  x  compo¬ 
nents  of  the  propagation  vectors  in  the  incident  beam 
and  Fix)  is  a  lateral  form  factor  determined  by  the 
transverse  shape  and  width  of  the  beam. 

By  Fourier  decomposition  the  field  B 3(x)  can  be  ex¬ 
pressed  as  an  integral  over  its  plane  wave  components 
associated  with  distinct  values  of  0  -  (c/u)k  j  in  the 
standard  form 

B3(x)  *  J*  d(363(/3)  exp  |i  —  0xj  (30) 

with 

b3(B)  “  y- J'  dxF(x)  exp  |i(fcf  -  —  fij  xj  •  (31) 

Each  plane  wave  component  in  B3(x)  contributes  a 
plane  wave  term  to  the  leflected  field  B3(x)  in  medium 
3  and  to  the  transmitted  field  Bj(x  ,z)  in  medium  1  de¬ 
termined  by  the  respective  plane  wave  reflection  and 
transmission  coefficients  r3210(/ 3)  and  a((3^)tz 21(d),  and 
the  total  fields  B3(x)  and  B\(x,z)  can,  therefore,  be  ex¬ 
pressed  in  the  forms 


which  is  easily  satisfied  in  the  cases  considered 
below. 

Given  that  relation  (37)  is  satisfied,  the  functions 
and  at321  in  the  integrands  of  Eqs.  (32)  and  (33)  [and. 
therefore,  (34b)  and  (35b)J  can  be  well  approximated 
by  their  single-pole  forms  defined  by  Eqs.  (21)  and  (27). 
and  (with  Imdm)  constrained  by  causality  to  be  posi¬ 
tive)  the  Fourier  transforms  R.mo  and  T32 1  can  be  re¬ 
duced  to  the  expressions 

2tt 

rt32lo(*  *“  *')  *  rwi&n) - 6(x  -  x') 

Ui/C 

1381 

T32 1<X  -  x',z)  -  tmidmj){m(x  -x').  (39) 

where 

Ui  V 

2irexp  i  —  om(x— x')%  x'  <  x. 
f«.(x  -  X')  a  •  C  (40) 

0,  x'  >  x. 

In  this  case  by  use  of  Eqs.  (34a)  and  (35a)  the  fields 
B;t(x)  and  Bjfx.o)  assume  the  forms 


B3(x)  * 

4. 

d3r32lo(B)b3(0)  exp  |i  j  |3xj  > 

(32) 

B  ,lx )  -  r:VJ<dm  1  +  2  -  >in2o  .im  1 

c 

Bi(x.y)  ■ 

f  dda(3-rlf32i(3)63(3)exp|i-3xj  • 

(33) 

X  J'  dx’B3 lx  )  exp  i  — 3»,<x-x') 

1411 
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BiU.2)  =  r  dx'Ba(i')  exp  i-^m(x-x')  • 

c  c 

(42) 

with 

.  ,  .  1  ,  1 

si n2orJJm>  =  -  rv,tdm) - - —  • 

■2  ry,(3n) 

The  first  term  on  the  right  in  Eq.  (41)  represents  the 
field  produced  by  the  direct  reflection  of  the  incident 
field  B;j  at  che  3-2  boundary,  which  would  occur  in  the 
absence  of  the  additional  2-1  and  1-0  boundaries,  while 
the  second  term  on  the  right  accounts  for  the  influence 
of  the  latter  boundaries  and  can  be  identified  as  the  field 
which  is  reradiated  back  out  of  the  boundary  layer 
mode. 

The  rapidly  varying  x  dependences  of  the  fields  B3, 
B and  Bj  derive  from  rapidly  varying  x -dependent 
phases,  which  are  approximately  determined  by  the 
central  k -vector  k?of  the  incident  field  B3.  As  a  con¬ 
sequence,  by  extracting  the  rapidly  varying  phase  factor 
expfik'lx)  from  the  fields,  a  set  of  more  slowly  varying 
field  functions  can  be  obtained  which  have  x  depen¬ 
dences  that  more  closely  characterize  the  more  impor¬ 
tant  spatial  variations  of  the  fields.  In  what  follows  we 
denote  the  more  slowly  varying  field  functions  by  a  tilde 
symbol  and  connect  them  to  the  fields  of  Eqs.  (41)  and 
(42)  by  the  relations 

B3(x)  =  S3(x)  exp<-ik°x),  S3(x)  =  B3(x)  exp(-ik9x), 

Bi(xu)  =  Bi(xx)  exp(-ikfa).  (43) 

Following  Ref.  3  the  integral  Eqs.  (41)  and  (42)  can 
be  replaced  by  coupled  (differential)  equations  for  the 
fields  B3  and  Bx  in  terms  of  the  incident  field  B3.  In 
particular,  by  differentiation  of  Eq.  (42)  with  respect 
to  x  and  use  of  Eq.  (42)  in  (41),  we  derive  two  coupled 
equations,  which  in  terms  of  the  slowly  varying  field 
functions  Bt(x,z)  and  B3(x)  assume  the  forms 

^  a, +-?*&,  <44) 

ax  c  c 

S3  *  r32^3  +  2  -r — sin(2032)fii,  (45) 

£  m 

where  2  (c/a>)k? and  r:)2,  tm,  and  <Ps2  are  evaluated  at 
d  =  dm- 

It  is  convenient  to  introduce  the  notation 

=  3'  +  iJ' 

«  =  — I  "  <3m  “  --(13"  4-  1  mO/lm) 

C  c 

ti-IR e3„  ~*0>,  146) 

( 

-  ^  '  J  -  - 

~ - Sin^032, 

C 

where  (w/elB'  and  (u7c)B"  represent,  respectively,  the 
real  and  imaginary  parts  of  the  propagation  vector  of 
the  mth  boundary  layer  mode  in  the  absence  of  the  3-2 
boundary  (with  the  mode  label  m  suppressed).  In  this 
notation  Eqs.  (44)  and  (45)  assume  the  simplified 
forms 


(44) ' 

(45) ' 


In  regions  of  x  beyond  the  point  xq  where  the  incident 
field  at  the  3-2  surface  terminates,  B3(x)  *  0,  and  the 
solutions  B3  and  of  Eqs.  (44)'  and  (45)'  have  x  de¬ 
pendences  expressed  by  the  simple  exponential 
exp[-a(x  -  xo)].  The  latter  expresses  the  exponential 
decay  of  the  boundary  layer  mode  whose  radiation  ac¬ 
counts  for  the  total  field  B3  at  values  oi  x  >  xq.  The 
decay  of  the  boundary  layer  mode  reflects  energy  losses 
from  the  mode  resulting  from  both  energy  radiation  into 
medium  3  and  energy  dissipation  in  medium  1.  The 
separate  losses  contribute  the  separate  positive  terms 
(w/c)  Imodm  and  (w/c)  Imdm  to  (w/c)  Im)3-n.  which  ap¬ 
pears  as  the  positive  real  part  of  a.  In  the  long-range 
surface-plasmon  mode  emphasized  in  the  present  paper, 
the  energy  losses  due  to  dissipation  are  minimized. 


1.  Case  of  Square  Wave  Input 

In  the  present  section  we  evaluate  the  fields  B3  and 
B 1  in  the  example  case  of  an  incident  beam  in  medium 
3  which  has  a  profile  as  a  function  of  x  that  closely  ap¬ 
proximates  a  square  wave.  In  this  case  the  lateral  form 
factor  F(x)  in  Eq.  (29)  can  be  represented  as  the  dis¬ 
continuous  function 


u)  w 

A. - <  x  <  - 

2  2 


where  w  defines  the  length  of  intersection  of  the  beam 
along  the  x  axis  [as  shown  in  Fig.  1(b)],  Substitution 
of  (29)  and  (47)  into  Eqs.  (41)  and  (42)  and  evaluation 
of  the  required  integrals  lead  to  reflected  and  mode 
fields  B3  and  Bi  of  the  forms 


n 

S3(x)  =*  r3263(x) - -  A  I 

a 


where  tv  and  tVi  are  defined  in  Eqs.  (46). 

It  follows  from  Eqs.  (49)  and  (46)  that  the  magnitude 
of  the  field  in  medium  1  is  maximized  for  a  given 
boundary  layer  mode  m  where  the  incident  beam  is 
tuned  so  that  the  central  value  of  It  ,  is  equal  to  ,w/c  ) 
Re^.  In  this  phase-matching  condition  (as  pointed 
out  in  Ref.  3),  the  oscillating  phase  factor  in  the  inte¬ 
grand  of  Eqs.  (42)  and  (41 )  is  nearly  stationary,  and  the 
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coupling  between  the  input  field  and  the  interface  mode 
is  optimized. 

For  a  tuned  input,  the  ratio  (2ai)/a  is  predominantly 
real  and  positive,  which  means  that  for  values  of  x  be¬ 
tween  -  w/2  and  w! 2  the  two  terms  in  Eq.  (48)  are  of 
opposite  sign.  The  resulting  cancellation  between  the 
two  terms,  in  certain  conditions  can  produce  an  inter¬ 
ference  zero  in  the  reflected  field  in  the  region  of  x  where 
the  incident  beam  intersects  the  3-2  boundary.  In 
these  conditions,  measurement  of  the  position  of  the 
zero  along  with  analysis  of  the  profile  of  the  reflected 
beam  provides  information  about  the  excited  interface 
mode,  which  cannot  be  extracted  from  measurements 
of  the  total  "reflectivity.  In  particular  by  measurement 
of  both  the  position  of  the  zero  in  the  reflected  profile 
and  the  rate  of  exponential  decay  of  the  tail  of  the  pro¬ 
file  (in  the  region  x  >  w/2),  the  independent  loss  terms 
in  the  mode  propagation  constant  associated  with  en¬ 
ergy  reradiation  and  dissipation  can  be  determined 
separately.  (In  contrast  a  measurement  of  the  width 
of  the  Lorentzian  dip  in  the  reflectivity  as  a  function  of 
H,  provides  only  a  determination  of  the  total  loss  term 
Re«.) 

Details  of  the  analysis  required  to  determine  the 
separate  loss  constants  lm5dm  and  lmdm  are  based  on 
Eq.  (48),  which,  for  x  in  the  interval  -w/2  <  x  <  w/2, 
gives  the  result11 


6'^x> 


1  -  exp 


(51) 


For  the  right-hand  side  of  Eq.  (511  to  vanish  at  a  value 
of  x  =  x  i  in  the  given  interval,  it  is  necessary  that 


~ -nhi)  »  rre,  (52) 

a 

where 

n(x)  s  ji  -  exp  (-«(*  +  7)])  ' 

With  the  beam  tuned  to  a  resonance  zero  in  the  reflec¬ 
tivity,  the  quantity  a  is  purely  real  and  equal  to  (w/c) 
lm(3m,  which  value  can  be  determined  by  measurement 
of  either  the  width  of  the  Lorentzian  dip  in  the  reflec¬ 
tivity  or  the  decay  rate  of  the  exponential  tail  in  the 
reflected  profile  (for  x  >xo  =  w/2).  Using  Eqs.  (46),  the 
real  and  imaginary  parts  of  Eq.  (52)  give  the  two  equa¬ 
tions 


—  Im6dm  - lm 

c  2n(xi) 

u>  or 

-  Re6dm  - Re 

c  2n(xi) 


f  32<flm  )  1 

sin2<i>i2<(L,  )| 
[  r32(i3m ) 
Isin2i»n2<d^,) 


(53a) 

(53b) 


where  n(xi)  can  be  computed  from  the  measured  value 
of  *[.  The  remaining  factor  on  the  right-hand  side  of 
Eqs.  (53)  can  be  rewritten  in  the  form 


1 


rntfm  > 


-<(1  -  l/r&(A, 


2  lsin2<tat(|3«  >1 
with  r:j2(/?m )  defined  by  the  explicit  formula 


<:\  ( 


7^ 


<  IV  <2  ' 


<jV  7|  -  3'm  +  <:iv  <1  -  3'i 


(54) 


(55) 


Since  (■>  and  <:t  are  presumed  known,  an  evaluation  of 
Eq.  (55)  requires  only  the  values  of  the  real  and  imagi¬ 


nary  parts  of  3m.  But  the  pole  approximation  under¬ 
lying  the  present  calculation  assumes  the  slowly  varying 
character  of  r32,  on  the  basis  of  which  r3 2<dm)  must  be 
approximately  equal  to  r32($m).  Therefore,  to  the  ac¬ 
curacy  of  the  present  calculation,  3m  can  be  replaced 
on  the  right-hand  side  of  Eqs.  (53)  and  (54)  by  3m,  the 
imaginary  and  real  parts  of  which  are  determined  by  the 
measured  value  of  a  and  the  tuning  angle  of  the  mode, 
respectively: 

3m  "  Red*  +  1  lm<3„  -  vTj  un*i  +  1  —  u  (56) 

Numerical  evaluation  of  Eqs.  (53)  then  allows  evalua¬ 
tion  of  the  loss  constant  lm6dm  and  pole-shift  R*6d«, 
associated  with  the  mode  radiation. 


III.  Numerical  Results  and  Analysis 

In  Figs.  2  and  3  we  plot  the  normalized  reflected 
power  Pr/P,  vs  distance  x  along  the  3-2  boundary  for 
the  case  of  a  square  wave  He-Ne  beam  incident  from  a 
prism  onto  a  dielectric-silver-dielectric  interface  with 
the  central  wave  vector  of  the  beam  tuned  to  the  long- 
range  surface  plasmon  mode.  The  plotted  curves  are 
derived  from  the  formula 


P. 

K  "  |flj(0)|*  ' 


(57  > 


1 

5 

C 

* 


*091 tt«* 

Fig.  2.  Ratio  of  reflected  power  to  incident  peak  power  plotted  aa 
function  of  transverse  distance  along  interface  measured  in  microns 
Position  zero  marks  center  of  incident  beam  which  extends  from  -500 
to  +.500  pm.  Parameters  are  as  listed  in  text  (a)  <fj  »  1.25  sm.  L/W 
-  0.35;  tb)  dj  -  1 11m.  L/W  -  0.23;  (c)  d:  -  0.9  am.  L/W  -  0.17*. 


Position 

Same  as  Fig.  2.  (d)  d2  »  0.75  (in.  L/W  «  0.106;  (e)  d3  •  0.6 
Mm.  L/W  »  0.059. 
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with  the  ratio  fijlx  )/&  »t0>  computed  from  Eqs.  U8)  and 
150)  and  the  vacuum  wavelength  A  and  dielectric  con¬ 
stants  taken  to  be  A  *  0028  A,  f„  ■  t:  »  1 15)-'.  « s  ■ 
11.72)*,  t[  *  <a«  *  -15.996  +  0.5256.  For  the  chosen 
parameters  and  an  assumed  metal  film  thickness  dt  of 
200  A  the  required  angle  of  incidence  is  ~62#. 

The  distinct  curves  in  Figs.  2  and  3  correspond  to 
different  values  of  the  cover  layer  thickness  d.  with  the 
beamwidth  u  fixed  at  l  mm.  The  curves  graph  the 
distorted  power  profiles  produced  by  the  surface  plas- 
mon  mode,  in  the  absence  of  which  the  reflected  power 
profile  must  closely  match  the  square  profile  of  the  in 
cident  beam.  The  distorted  profiles  result  from  a 
•  partial)  absorption  of  the  energy  of  the  incident  beam 
into  the  mode  and  the  subsequent  reradiation  of  this 
energy  back  into  the  incident  medium.  Beyond  the 
edge  of  the  region  illuminated  by  the  incident  beam  (for 
x  >  xo  ■  w/2),  where  the  reflected  power  P,  derives 
strictly  from  the  reradiated  plasmon  field.  P,  decreases 
exponentially  with  a  decay  constant  a  given  by  (w/ci 
lm3,*  On  the  other  hand,  for  values  of  x  within  the 
illuminated  region  of  the  surface  (for -w/2  <  x  <  w/2). 
where  the  reflected  power  results  from  both  the  rera- 
diatied  field  and  the  directly  reflected  incident  field, 
destructive  interference  between  these  two  fields  results 
in  a  distorted  square  profile,  details  of  which  are  de¬ 
termined  for  given  dielectric  parameters  by  the  propa¬ 
gation  length  of  the  surface-plasmon  mode  in  relation 
to  the  width  of  the  incident  beam.  The  shift  of  the 
centroid  of  the  reflected  profiles  in  Fig.  3  relative  to  that 
of  the  incident  profile  corresponds  to  the  so-called  Goos 
Hanchen  shift.812 

The  propagation  length  L  of  the  surface-plasmon 
mode,  defined  here  as 


L  » 


l 

<*i|c  Imi?, 


“  di'  +  Imiid^) 

c 


(58) 


in  the  reflected  profile  in  the  region  illuminated  by  the 
beam  To  obtain  a  profile  with  a  well  defined  zero  »uch 
as  shown  in  curve-  u i  (  and  iooI  Fig  3  the  ratio  of  the 
propagation  length  L  to  the  beamwidth  a  needs  to  be 
restricted  between  definite  limits  of  less  than  one  In 
the  particular  case  shown  in  Figs  2  and  3  the  ratio  L  u 
needs  to  be  restricted  between  the  limns  1  5  and  1  2() 
For  larger  value*  of /.  u  the  intensity  of  the  reradtated 
plasmon  field  is  t.*.  weak  in  the  region  illuminated  by 
the  beam  to  produce  an  interference  zero  in  that  region, 
while  for  smaller  value*  of  L  u  the  interference  zero  can 
occur  too  close  to  the  leading  edge  of  the  reflected  beam 
to  be  observable  In  practice,  for  a  given  input  beam 
width,  given  ihat  the  plasmon  low*  in  the  metal  are 
sutfcientlv  small,  an  interference  zero  in  the  reflected 
profile  can  he  obtained  bv  adjustment  of  the  cover 
thickness.  In  this  case  analysis  of  the  profile  described 
in  Sec.  II  allows  for  determination  of  the  separate  ’.os* 
coefficients  d*  and  lmo.i 


Appendix 

In  this  Appendix  we  give  an  expression  for  the 
quantity  (dv  >/<dd)  in  Eq.  <  181  and  present  a  general 
formula  for  the  power  per  unit  length  P  carried  by  an 
interface  mode.  The  result  for  P  is  valid  for  complex 
dielectric  constants  in  media  2.  1.  and  0  and  represents 
a  generalization  of  the  formula  given  by  Marcuse  in 
Ref.  1. 

We  consider  the  case  of  mode  excitation  by  TM 
polarized  radiation  with  longitudinal  k -vector  k  * 
(w/c)i3.  In  this  case,  with  <b,i  *  tan_'U;fc|X/uifc„ )  and 
i£(/5)  *  2ikitd\  -  O12  -  <J>iob  we  find  for  the  quantity 
(d\£)/(di3)  the  result 


where 


d3 


'AD 


is  determined  by  the  (fixed)  metal  film  thickness  di, 
which  determines  the  loss  coefficient  due  to  dissipation 
0"  and  by  the  (varied)  cover  thickness  do,  which  de¬ 
termines  the  loss  coefficient  due  to  reradiation  lm5|5. 
Where  the  propagation  length  of  the  mode  is  short 
compared  with  the  width  of  the  beam,  the  plasmon  field 
excited  by  the  leading  edge  of  the  beam  is  reradiated 
strictly  into  the  region  illuminated  by  the  beam.  In  this 
case  the  destructive  interference  between  the  reradiated 
field  and  the  directly  reflected  field  can  produce  a  zero 
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An  interface  mode  of  the  (free)  layered  structure 
consisting  of  the  2-1  and  1-0  boundaries  consists  in 
general  of  multiply  reflected  fields  in  medium  1  with 
/r-vectors  k,  and  k|  and  outgoing  fields  in  media  2  and 
0-with  k -vectors  kj  and  k,(  as  sketched  in  Fig.  4.  The 
time  averaged  power  flow  carried  in  the  mode  along  x 
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per  unit  length  in  they  direction  is  computed  from  the 
formula 

P=£  dz(S)  x  =  -^-Re  j"  dzE,B’y,  (A3) 

where  (S)  represents  the  time  averaged  Poynting  vec¬ 
tor.  The  power  P  is,  therefore,  the  sum  of  the  power  per 
unit  length  associated  with  the  mode  fields  in  the  three 
regions  0,  1,  and  2  (in  the  absence  of  a  3-2  boundary). 

By  computing  the  mode  fields  in  the  three  media 
subject  to  the  boundary  matching  relations  and  the 
mode  condition  </>i2  +  <2>io  =  &irdi  -  mir  and  performing 
the  required  integration,  we  obtain  for  P  the  generad 
result 


In  the  case  of  either  a  dissipationless  surface-plasmon 
mode  where  lme;  =  0,  lm/3  =  0,  R ekjt  *  0,  for  j  *  0,1,2 
or  a  dissipationless  guided  wave  mode  where  (for  j  - 
0,1,2)  lme;  =  0,  lmd  =  0,  R efe0*  =  Refe2 z  ~  0,  and  1mA  u 
=  0,  Eq.  (A4)  for  P  reduces  to 
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surface  plasmon, 
guided  wave. 
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Xhe  result  ( A7)  is  consistent  with  the  result  given  in  Ref. 

1. 


P  =•  —  |a 1 1 2( RetJ)  exp<-2  Im0,2  -  2  lmA^xli  |, 
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where  a  t  is  the  amplitude  of  the  plane  wave  field  with 
£ -vector  k[  in  medium  1,  and  the  bracket  is  expressible 
in  either  of  the  forms 
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Optics  ’83:  A  Report  on  Emerging  Technologies 

The  2nd  Annual  Report  on  Optics  by  the  Optical  Society  of  America's 
Technical  Council  profiles  research  findings  from  around  the  country. 


Long-Range  Surface- 
Plasmon  Polaritons 

The  physics  of  surface-plastnon 
polaritons  propagating  on  very  thin  metal 
films  that  are  bounded  by  similar  and 
semi-infinite  dielectrics  or  semiconduc¬ 
tors  has  been  extensively  studied  in  the 
past  two  years.  The  interest  in  these 
polaritons  arose  after  it  was  discovered1-2 
that  they  are  accompanied  by  large  elec¬ 
tric  and  magnetic  fields  that  can  serve  as  a 
vehicle  for  enhancing  linear  and  nonlin¬ 
ear  interactions.  These  polaritons  consist 
of  two  modes  that  can  be  excited  prefer¬ 
entially.  One  mode  has  a  range  along  the 
metal  film  that  increases  as  the  metal  film 
thickness  on  which  it  propagates  de¬ 
creases  (Long-Range  Surface-Plasmon, 
or  SRSP).  The  other  mode  has  the  oppo¬ 
site  property,  namely,  that  this  range  de¬ 
creases  as  the  metal  film  thickness  de¬ 
creases  (Short-Range  Surface-Plasmon, 
.gw  or  LRSP).  The  range  of  a  mode  is  deter- 
mined  by  the  ratio  of  the  power  that  is 
carried  by  the  lossy  metal  film  and  by  the 

lossless  bounding  media  via  the  evanes¬ 
cent  fields:  For  the  LRSP  (SRSP),  as  the 
metal  film  thickness  decreases,  the  frac¬ 
tion  of  the  power  carried  by  the  evanes¬ 
cent  fields  increases  (decreases),  resulting 
in  an  increase  (decrease)  in  the  range.* 

The  first  experimental  observation  of 
the  difference  in  the  properties  of  these 
two  modes  was  reported  by  Kovacs1.  The 
first  experimental  tracing 4  of  the  theoret¬ 
ical  dispersion  curves',  and  the  fit  be¬ 
tween  the  calculated  and  measured  Lorent- 
zian  absorption  curves1  of  these  two 
modes  were  recently  reported  with  film 
thicknesses  ranging  from  100  nm  down  to 
10  nm. 

There  are  several  important  interac¬ 
tions  that  can  take  advantage  of  the  en¬ 
hanced  fields  of  the  LRSP.  For  example, 
a  recent  calculation  involving  a  single 
interface  surface-plasmon  polariton  has 
revealed  that  it  can  have  a  propagation 
constant  which  depends  on  the  optical 
power  that  it  carries4.  By  considering  the 
two-interface  polaritons,  these  two  the¬ 
ories1-6  could  now  be  merged.  It  was 
found  that  the  LRSP  indeed  has  an  en¬ 
hanced  optical  power  dependence  of  its 
propagation  constant7-8.  We  are  currently 
working  on  experiments  to  observe  bista¬ 
bility  in  a  system  which  utilizes  the  LRSP 
propagating  on  thin  metal  films  bounded 
by  two  InSb  crystals.’ 


Sacond-harmonic  generation 
Another  application  of  the  enhance¬ 
ment  in  the  optical  fields  that  accompany 
the  LRSP  is  second-harmonic  generation. 
It  has  been  calculated  and  observed  ex¬ 
perimentally  that  single-interface  polar¬ 
itons  can  give  rise  to  an  enhancement  in 
second  harmonic  generation  on  total  in¬ 
ternal  reflection10.  One  could  now  use 
two-interface  polaritons  and  show  the¬ 
oretically  that  the  LRSP  will  give  an  in¬ 
crease  of  more  than  two  orders  of  magni¬ 
tude  in  this  enhancement' ' .  The  predicted 
increase  in  the  enhancement  of  the  second 
harmonic  generation  efficiency  has  been 
recently  observed  experimentally'  *. 

We  expect  that  the  LRSP  will  be  found 
useful  for  other  applications  as  well.  Be¬ 
cause  of  the  ease  with  which  thin  metal 
films  can  be  deposited  on  practically  any 
substrate,  these  modes  can  be  used  to 
probe  metal-dielectric  interface  properties 
such  as  roughness,  or  charge  carrier  prop¬ 
erties  in  very  thin  metal  films,  for 
example.  -  Dror  Sarid,  Optical  Sciences 
Center,  University  of  Arizona. 
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ABSTRACT 

A  novel  system  for  coupling  to  surface  plasmon  polaritons  has  been  designed, 
fabricated  and  assembled.  Microcomputer  controlled  electronics  provide  for  accurate 
and  repeatable  angular  positioning  of  the  waveguide  coupler  with  respect  to  the  incident 
beam,  and  for  precise  angular  registration  of  detector  output  signals.  The  optical 
components,  combined  with  a  retroreflecting  spherical-surface  coupling  prism,  give 
apertured  control  of  the  beam  profile  at  the  coupling  interface. 
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INTRODUCTION 


Prism  couplers  have  been  used  extensively  to  provide  evanescent  field  (frustrated 
total  internal  reflection)  coupling  to  optical  guided  wave  structures.  A  discussion  of 
waveguide  couplers  has  been  presented  by  P.  K.  Tien1.  Prism  coupler  systems  usually 
comprise  a  triangular  cross- sectioned  prism  of  high  refractive  index  glass  mounted  on  a 
rotary  table  secured  by  one  of  a  variety  of  clamping  structures  which  also  controls  the 
proximity  of  the  prism  to  the  waveguide.  Properly  constructed  clamping  schemes  can 
produce  coupling  efficiency  as  high  as  94%  2,3.  The  entire  assembly  is  positioned  so 
that  the  center  of  rotation  of  the  table  and  the  center  of  the  prism-waveguide 
combination  on  the  table  lie  in  the  path  of  a  quasi-collimated  beam  of  a  laser.  The 
table  is  rotated  until  the  condition  for  coupling  the  incident  light  to  the  waveguide 
structure  is  satisfied. 

To  observe  either  the  coupling  efficiency  or  the  angular  width  of  the  coupling 
regime,  the  power  of  the  beam  which  is  reflected  from  the  coupling  interface  is 
monitored  by  a  detector  as  a  function  of  the  angle  of  incidence.  For  a  symmetrically 
shaped  triangular  prism,  the  reflected  beam  exits  the  prism  at  an  angle  relative  to  the 
coupling  interface  which  is  equal  to  the  angle  between  the  incident  beam  and  the 
interface.  Thus  the  angle  between  the  incident  and  reflected  beams  is  twice  the  angle 
separating  the  incident  beam  from  the  normal  to  the  coupling  interface.  Rotating  the 
prism  causes  the  reflected  beam  direction  to  change  at  an  angular  rate  twice  that  of  the 
prism.  A  description  of  the  propagation  direction  of  the  reflected  beam  for  an 
asymmetric  triangular  prism  requires  a  more  complicated  calculation.  Arranging  for  this 
light  to  fall  continuously  on  a  detector  as  the  table  and  prism  waveguide  structure 
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rotate  i$  difficult  with  systems  incorporating  triangular  prisms. 

Triangular  prism  systems  possess  other  undesirable  optical  characteristics.  As  the 
prism  is  rotated,  the  point  where  the  incoming  beam  intersects  the  coupling  interface 
walks  along  the  interface.  There  is  no  axis  around  which  the  prism  can  be  rotated  to 
ensure  that  the  beam  always  intersects  the  coupling  interface  at  a  fixed  point.  Unless 
the  low  refractive  index  gap  region  across  which  the  light  couples  to  the  waveguide  is 
perfectly  uniform  in  thickness,  this  walking  causes  variation  in  the  coupling  efficiency. 
The  efficiency  is  especially  sensitive  to  gap  thickness  for  coupling  angles  near  the 
critical  angle  in  the  prism.  Moreover,  loss  by  reflection  as  the  light  enters  the  prism 
varies  with  angle  of  incidence,  according  to  the  Fresnel  equations.  These  factors  cause 
variations  in  the  measured  signal  that  are  not  characteristic  of  the  guided  mode,  and  for 
which  compensation  can  be  difficult  or  even  impossible. 

The  system  described  below  resolves  many  of  the  undesirable  features  of  traditional 
prism  coupler  systems.  In  addition,  it  provides  flexibility  for  controlling  the  beam's 
amplitude  profile  and  transverse  extent  at  its  intersection  with  the  coupling  interface, 
which  is  particularly  important  when  measuring  the  shape  of  the  profile  of  the  reflected 
beam.4 


OPTICAL  SYSTEM 

A  schematic  of  the  optical  system  is  shown  in  Fig.  1,  and  a  photograph  and  sketch  of 
the  prism  appear  in  Fig.  2.  This  prism  is  the  keystone  of  the  design. 

The  prism  is  made  of  a  high  index  glass,  SF56,  having  a  refractive  index  at  a 
wavelength  of  632.8  nm  of  n=1.7788  (interpolated  from  Schott  catalog  data).  Before 
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grinding  the  hemispherical  surface,  the  glass  block  was  polished  flat  on  one  side,  sawed 
in  half,  mirror  coated  with  aluminum  on  the  polished  sides,  then  cemented  together  with 
the  mirrored  sides  in  contact.  Next,  a  base  was  polished  perpendicular  to  the  internal 
mirrored  surface  forming  a  retroreflecting  corner.  Then  the  hemispherical  surface, 
centered  on  the  intersection  of  the  base  with  the  internal  mirrored  surface  was  polished. 
The  final  radius  of  the  hemisphere  is  11.8  mm. 

The  micrometer  pressing  on  the  base  of  the  prism  controls  the  thickness  of  the  low 
index  gap  across  which  the  light  couples  to  the  waveguide.  By  adjusting  this 
micrometer,  the  amount  of  light  coupled  to  the  mode  can  be  varied.  The  coupling 
constant  is  given  by4 

h  =  exp^-(k,-^-)1/*  d,j,  (1) 

where  k  is  the  longitudinal  component  of  the  wave  vector,  c  is  the  speed  of  light  in 
vacuum,  u>  is  the  frequency  of  the  light,  e2  is  the  dielectric  constant  of  the  gap  and  d2  is 
the  thickness  of  the  gap.  In  our  system  the  coupling  constant  can  be  varied  from  0  to 
1/2. 

The  beam  from  the  HeNe  laser  is  incident  on  a  spatial  filter  (a  3.8  mm  focal  length, 
infinite  conjugate  microscope  objective  with  a  5  micron  pinhole  at  its  focal  point),  and 
an  achromatic  doublet  lens  collimates  the  light.  An  adjustable  rectangular  aperture  and 
a  beam  splitter  are  located  in  the  collimated  space.  A  second  achromatic  doublet  brings 
the  beam  to  a  focus  at  the  front  focal  point  of  the  spherically  surfaced  prism,  whose 
base  is  the  coupling  interface.  The  spherical  surface  again  collimates  the  light.  The  ray 
bundle  will  retroreflect,  bouncing  first  off  the  internal,  mirrored  surface,  then  off  the 
coupling  surface  (or  perhaps  in  reverse  order).  These  rays  then  traverse  the  hemisphere 
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again  and  are  refracted  to  pass  through  the  same  focal  point.  Now  ro  ating  the  prism 
about  its  polar  axis  (defined  by  the  retroreflecting  corner)  causes  the  angle  of  incidence 
of  the  light  bundle  on  the  coupling  interface  to  change  by  the  same  amount.  Since  the 
retroreflecting  corner  and  the  center  of  the  hemisphere  lie  on  this  axis,  the  beam  does 
not  walk  along  the  coupling  interface  except  to  the  extent  that  its  intersection  width 
increases  as  a  result  of  its  incident  obliquity.  In  addition,  the  reflected  beam  is  always 
found  in  the  same  place  and  propagating  in  the  same  direction.  Although  curved  surface 
prism  couplers  have  been  used  by  other  researchers^  this  one  is  unique  in  that  it  utilizes 
a  retroreflecting  hemispherical  design. 

The  rectangular  aperture  is  located  such  that  it  is  imaged  by  the  second  achromat 
and  the  spherical  surface  of  the  prism  onto  the  coupling  interface.  The  demagnification 
of  these  two  elements  is  the  ratio  of  their  focal  lengths,  1:2.65.  In  use,  the  aperture 
dimensions  are  generally  about  4  mm  square;  the  imaged  dimensions  are  thus  about  1.5 
mm.  Adjusting  the  aperture  allows  different  portions  of  the  beam  to  pass  for  sampling 
various  regions  of  the  coupling  interface.  Furthermore,  the  beam  incident  on  the 
interface  maintains  its  rectangular  intensity  profile.  Care  must  be  taken  that  the  slit 
dimensions  not  be  made  too  small,  since  the  demagnification  increases  the  divergence  of 
the  beam's  angular  spectrum. 

The  beam  splitters  direct  half  of  the  return  beam  from  the  prism  coupler  to  two 
detectors.  The  first  beam  splitter  also  directs  some  of  the  light  to  a  slit  and  a  mirror 
used  for  angle  calibration.  Light  passes  through  the  slit  to  the  mirror  and  off  the 
substrate  holding  the  waveguide.  At  a  well  defined  angle,  the  light  will  return  through 
the  slit  to  reach  the  detectors.  This  angle  then  provides  a  reference  from  which  other 

•  1 

directions  can  be  measured  accurately.  The  calibration  is  reproducible  to  degree. 
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One  might  expect  the  spherical  surface  of  the  prism  combined  with  the  prism's  high 
refractive  index  to  introduce  a  severe  amount  of  positive  spherical  aberration  to  the 
beam.  This  is  not  the  case  in  our  system,  since  we  use  a  small  aperture.  Standard  texts 
treat  aberrations  in  detail.  In  our  case,  the  spherical  aberration  is  given  by 


W 


1  y*n'*(n'-1 ) 


«m  a  r 


p -  (2) 

where  y  is  the  marginal  ray  height  at  the  surface,  n1  is  the  index  of  refraction  of  the 
prism,  and  r  is  the  radius  of  curvature  of  the  surface.  With  y=i  mm,  r=11.8  mm  and 
n'si.7788  a  simple  calculation  shows  that  the  aberration  introduced  is  =  0.19  um 
(less  than  a  third  of  a  wavelength). 

If  this  aberration  is  a  problem,  additional  correcting  elements  may  be  added  to  the 
system.  In  particular,  a  negative  meniscus  lens  may  be  inserted  between  the  spatial 
filter  and  the  collimating  lens  to  introduce  negative  spherical  aberration  to  the  beam, 
exactly  compensating  for  the  aberration  of  the  prism.  The  collimating  lens  must  be 
repositioned  to  maintain  good  collimation,  and  adjusting  the  axial  position  of  the 
meniscus  lens  changes  the  amount  of  aberration  introduced.  The  aberration  correction 
can  then  be  checked  by  inserting  a  plane  parallel  plate  shearing  interferometer  into  the 
collimated  beam  to  simultaneously  measure  the  aberration  in  the  incoming  and  outgoing 
beams.  Since  the  prism  is  used  in  double  pass,  the  aberrations  in  the  incoming  and 
outgoing  beams  have  the  same  magnitude  (but  opposite  sign)  when  the  beam  inside  the 
prism  is  well  corrected.  The  meniscus  lens,  however,  introduces  additional  light 
scattering  from  flaws  and  dust  on  its  surface;  if  beam  uniformity  is  important  then  the 
interference  produced  by  this  scattered  light  may  cause  problems  more  severe  than  the 
spherical  aberration.  For  this  reason  the  aberration  control  was  removed  from  our 


system 
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ELECTRONIC  SYSTEM 

The  rotation  of  the  table  which  supports  the  hemispherical  waveguide  coupler,  and 
the  registration  and  recording  of  the  detected  light  signal,  are  controlled  via  a 
microcomputer  in  conjunction  with  various  accessory  electronic  instruments.  A  flow 
chart  of  the  control  chain  is  presented  in  Fig.  3. 

One  of  the  detectors  is  a  wavelength-filter-compensated  silicon  photojunction 
detector  operated  in  photoconductive  mode.  It  is  preceded  by  a  narrow  band 
wavelength  filter  which  allows  operation  when  the  room  lights  are  on  without  serious 
degradation  of  the  signal-to-noise  ratio.  The  detector  and  its  preamplifier  are  sensitive 
to  a  wide  range  of  incident  power,  from  0.01  yW  to  1  mW.  From  the  detector  the  signal 
passes  through  a  preamplifier  to  a  digital  multimeter.  The  signal  is  then  transmitted  in 
digital  form  to  the  microcomputer  by  the  DMM  when  it  is  queried. 

The  other  detector  is  a  CCD  1728  element  linear  image  sensor  positioned  so  that  the 
coupling  interface  is  imaged  onto  the  detector.  This  enables  one  to  examine  the 
intensity  profile  of  the  reflected  beam.  A  discussion  of  the  purpose  and  use  of  the  CCD 
array  will  be  presented  elsewhere. 

Completing  the  control  loop,  the  microcomputer  addresses  and  instructs  a  controller 
unit  which  produces  the  required  pulse  signals  to  drive  the  rotary  table,  on  which  the 
waveguide  coupler  is  mounted  and  centered,  to  the  desired  angular  position.  The  table 
is  driven  by  a  worm  gear  mechanism  and  may  be  positioned  repeatably  to  within  3  arc 
minutes  cumulative  accuracy  and  0.2  arc  minutes  repeatibility,  with  an  angular  step  size 

of  -1  arc  minute  (0.00417  degree). 
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APPLICATIONS  ANO  RESULTS 

Our  coupling  system  has  been  used  to  probe  surface  plasmon  polaritons  on  thin  silver 
films&»7,  using  the  attenuated  total  reflection  (ATR)  approach.  This  type  of  resonance  is 
associated  with  a  light  wave  guided  by  a  symmetrical  refractive  index  structure 
consisting  of  dielectric,  metal,  and  dielectric  layers.  The  attenuation  of  the  polariton 
along  the  metal  film  is  reasonably  low  (in  the  range  of  30  cm-1)  and  therefore  its 
resonance  is  fairly  narrow.  The  theoretical  width  of  the  resonance  is  0.001  degree  for  a 
silver  film  5  nm  thick  and  a  wavelength  of  632.8  nm.  For  the  tested  films,  as  thin  as  10 
nm,  the  angular  widths  of  the  waveguide  resonances  have  been  as  narrow  as  20%  of  the 
angular  resolution  achievable  by  the  optical  system,  allowing  close  qualitative  analysis  of 
the  system's  capabilities.  A  plot  of  the  angular  reflectance  of  a  very  narrow  resonance 
(theoretical  width,  0.012  degrees)  produced  by  the  system  is  shown  in  Fig.  4. 

The  rectangular  delimiting  of  the  beam's  cross-section  at  the  coupling  interface  is  a 
vital  part  of  the  system's  design.  It  facilitates  analysis  of  the  various  loss  components 
of  the  guided  mode,  which  consist  of  absorption  in  the  metal,  roughness  scattering  and 
out-coupiing  to  the  prism.4  However,  it  is  also  the  apparent  limiting  feature  of  the 
optical  system's  performance.  Diffraction  from  a  rectangular  aperture  only  1.0  mm  in 
extent  produces  a  beam,  at  wavelength  632.8  nm,  that  has  a  diffraction  angular 
spectrum  of  0.044  degrees,  FWHM.  An  observed  resonance  0.4  degree  wide  showed 
diffraction  ringing  when  the  aperture  width  was  reduced  below  0.8  mm.  (With  this 
aperture  width,  the  first  zero  of  the  sine  function  of  the  angular  spectrum  is  at  0.38 
degrees.)  This  angular  width  is  more  than  ten  times  the  resolution  of  our  rotary  table. 
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The  small  effective  aperture  is  the  result  of  the  demagnification  ratio  (1:2.65)  of  the 
last  achromat  and  the  hemispherical  coupler.  An  increase  in  the  optical  angular 
resolution  may  be  obtained  by  increasing  the  size  of  the  adjustable  rectangular  aperture. 
This  change  also  requires  greater  attention  to  control  of  the  off-axis  beam  quality. 
Decreasing  the  focal  length  of  the  second  achromat  would  also  help  but  is  not  practical 
because  of  the  extremely  short  focal  length  and  high  numerical  aperture  that  would  be 
required. 

Although  we  have  used  the  system  for  measuring  the  dispersion  curves  of  surface 
plasmon  polaritons  propagating  on  very  thin  metal  films,  the  system  can  be  used  for 
exciting  and  probing  optical  waveguides  and  other  surface  modes  as  well. 

CONCLUSION 


A  hemispherical,  retroreflecting,  ATR  waveguide  coupler  system  has  been 
implemented  v-nich  exhibits  several  attractive  features  compared  to  traditional  prism 
coupler  systems.  Its  use  in  a  microcomputer  controlled  system  to  measure  narrow 
angular  resonances  has  been  demonstrated. 
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FIGURE  CAPTIONS 

Optical  system  layout:  Support  optics  for  re tron? fleeting  hemispherical  coupler. 

a)  Photograph  of  the  coupler  with  retroreflectLng  laser  beam. 

b)  Sketch  of  an  overhead  view  of  the  coupler. 

Flow  chart  of  control  electronics. 

System  resolution  of  approximately  0.044  degree  is  apparent  from  this 
experimental  trace  of  a  narrow  (0.012  degree,  theoretically)  surface- plasmon 


ATR  resonance. 
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